THE SINE-GORDON EQUATION IN THE SEMICLASSICAL LIMIT: CRITICAL 

BEHAVIOR NEAR A SEPARATRIX 
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Abstract. We study the Cauchy problem for the sine-Gordon equation in the semiclassical Umit with pure- 
impulse initial data of sufficient strength to generate both high-frequency rotational motion near the peak 
of the impulse profile and also high-frequency librational motion in the tails. Subject to suitable conditions 
of a general nature, we analyze the fluxon condensate solution approximating the given initial data for small 
time near points where the initial data crosses the separatrix of the phase portrait of the simple pendulum. 
We show that the solution is locally constructed as a universal curvilinear grid of superluminal kinks and 
grazing collisions thereof, with the grid curves being determined from rational solutions of the Painleve-II 
system. 
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1. Introduction 



This paper is concerned with a detailed local analysis of the solution of the Cauchy initial-value problem 
for the sine-Gordon equation 

e^utt - e^Ua;x + sm{u) = 0, u{x,0) = F{x), eut{x,0) ^ G{x), xeM.. (1.1) 

We will consider the number e > to be a small parameter. This type of scaling can be physically motivated 
in the situation that the sine-Gordon equation is used to model the propagation of magnetic flux along 
superconducting Josephson junctions |23j . The sine-Gordon equation can also be derived in the continuum 
limit as a model for an array of coaxial pendula with nearest-neighbor torsion coupling [3]. This latter 
application is particularly useful for the purposes of visualization of solutions. 

A dramatic separation of scales occurs in the semiclassical limit e J, if the initial data F{-) and G{-) 
are held fixed. As can be seen in Figures |l.lf[1.3[ the semiclassical dynamics apparently consists of well- 
defined (asymptotically independent of e) spacetime regions containing oscillations on space and time scales 
proportional to e but modulated over longer scales originating with the e-independent initial conditions. An 
important role is played by the x-parametrized curve {F,G) = {F{x),G{x)) in the phase portrait of the 
simple pendulum, and its relation to the separatrix curve (1 — cos(F)) + ^G^ ~ 2. Indeed, in the specific 
context of suitable initial data of pure impulse type, that is, for which F{x) = 0, the following dichotomy 
has recently been established [7] regarding the asymptotic behavior of the solution u{x,t) of ( |1.1[ ) for small 
time t independent of e. If {F,G) = (0,G(x)) hes inside the separatrix, then u{x,t) is accurately modeled 
by a modulated train of superluminal librational waves; but if {F, G) = (0, G{x)) lies outside the separatrix, 
then u{x, t) is instead accurately modeled by a modulated train of superluminal rotational waves. If a; S M is 
a point lying exactly on the separatrix curve, then the approximation theorems proved in [7] fail to provide 
a uniform description of the asymptotic behavior of the solution u{x,t) near such x. It is clear from plots 
of exact solutions shown in Figures |1. 1^1.31 that some essentially different asymptotic behavior is generated 
by separatrix crossings in the pure-impulse initial data. In particular, a different and more complicated 
kind of waveform than modulated traveling waves appears to spread in time away from specific points where 
separatrix crossings occur in the initial data. The coupled pendulum interpretation is useful here: if the 
pendula are all initially at rest in the gravitationally stable configuration and are given a spatially-localized 
initial impulse of sufficient strength, then some pendula have sufficient energy to rotate completely about the 
axis a number of times, while the pendula in the "wings" experience very little initial impulse and only have 
energy for small oscillatory motions near equilibrium (so-called librational motion). Clearly this situation 
leads to kink generation near the transition points and a more complicated type of dynamics as the kinks 
struggle to separate from one another. 

To better understand the reason for the locally complicated behavior, it is useful to view the sine-Gordon 
equation as a perturbation of the simple pendulum ODE: 

UTT + sin(u) = e^u^x, u{x, 0) = 0, uxix, 0) = G{x), t = eT. (1.2) 

We might expect that the right-hand side should be negligible for quite a long rescaled time T as long 
as in the unperturbed system (e — 0) pendula located at nearby values of x follow nearby orbits of the 
pendulum phase portrait. This will be the case unless |G(a;)| « 2, the condition allowing for nearby values 
of X to correspond to topologically dissimilar orbits, leading to large relative displacements in u over finite 
T and causing Uxx to become very large very quickly. See Figure |1.4[ To leading order, and for times t of 
order e we expect that orbits near x — Xait should follow the pendulum separatrix: cos(u/2) « sech(i/e), 
sin(M/2) « — tanh(t/e), and eut ~ — 2sech(t/e). Of course as T = t/e — >■ oo, the perturbation term e^Uxx 



will become important and will prevent the rapid divergence of trajectories shown in Figure 1.4 



The subject of this paper is the asymptotic analysis, in the semiclassical limit e J, 0, of the solution of the 



Cauchy problem ( 1.1 ) for the sine- Gordon equation in the case that t is small and x is near a value where the 



curve {F,G) = {0,G{x)) crosses the separatrix (at its midpoint, as in Figure 1.4). This region of the {x,t)- 
plane is blown up with what turn out to be the correct scalings for easier viewing in Figure |1.5| Therefore 
we are near the boundary of the regions where the solution is described in terms of modulated rotational 
and librational solutions of the pendulum equation. We will show that, under suitable further assumptions, 
the asymptotic behavior of u{x,t) in this situation is universal, and is described by an essentially multiscale 
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Figure 1.1. Density plots of cos{u(x,t)) (left) and sm{u{x,t)) (right) for the exact solution 
of the Cauchy problem (1.1 ) with e = 0.1875 for F{x) = and G{x) = — 3sech(a;). For this 
data, separatrix crossings occur for x = ±arcsech(|) w ±0.962. 



formula that is expressible in terms of modern special functions, specifically solutions of certain nonlinear 
ordinary differential equations of Painleve type. 

To formulate our results precisely requires that we first set up some background material; we hope that 
the reader will bear with us until §1.3| where the full details will be presented. In the meantime, we can 
describe the semiclassical asymptotics of u{x,t) in the region illustrated in Figure 1.5 by saying that this 



region contains a curvilinear network of isolated kink-type solutions of the sine-Gordon equation (with the 
approximation u{x, t) « tt (mod 2tt) holding in between the kinks) lying along the graphs of certain rational 
solutions of inhomogeneous Painleve-II equations. At spacetime points associated with poles and zeros of 
these solutions, the kinks collide in a grazing fashion as can be modeled by a more complicated exact solution 
of the sine-Gordon equation corresponding to a double soliton. The phenomenon of an asymptotically 
universal wave pattern being described by simple "solitonic" solutions located in space and time according 
to a rather more transcendental solution of a Painleve-type equation was only recently observed for the first 
time in another problem by Bertola and Tovbis [5]. 

1.1. Notation and terminology. All power functions z^' will be assumed to be defined for nonintegral real 
p as the principal branch with branch cut z < and with — |p|7r < arg(zP) < \p\7r. 

If G(-) is a function of x for which the equation G{x) — —2 has a unique solution for x > 0, then we 
denote this solution by cc = Xcrit- Let us use the term criticality to describe the point x — Xcnt and t = 0, 
and near criticality to mean x ~ a;crit and t « 0. 

For any rational function TZ = Tl{y), let ^{TZ) denote the finite set of real poles of TZ and let 3f{TZ) 
denote the finite set of real zeros of TZ. Finally, for a finite set ^ C M and a real number y E M. denote by 

yo\ (1.3) 



\y 



.y \ :~ min \y - 



the distance between y and J^. 
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Figure 1.2. Same as Figure [7T7| hut with e = 0.09375. 



If S is an oriented contour in the complex plane and / is a function analytic in the complement of S, 
we will use subscripts f+{£,) and /-(^) to refer to the boundary values taken by f{w) as w — > ^ e S 
nontangentially from the left and right respectively. 

We will make frequent use of the Pauli matrices: 



0-1 



ro 


1' 




-i 




'1 " 


1 — ' 





, CT2 := 


i 


, 0-3 :== 


-1 



(1.4) 



With the sole exception of these three, we write all matrices with boldface capital letters. 

We use the Landau notation for most estimates, with "big-oh" written O and "little-oh" written o. Also, 
if 9i: 92, • ■ • , 9m are some quantities, we will use the shorthand notation 0{qi, . . . ,qm) to represent a quantity 
that is bounded by a linear combination of |(7i|, . . . , \ qm I, that is, 0{qi, . . . ,q,n) = 0{qi) + ■ ■ ■ + 0{q^). 

We will be dealing with several matrix functions in which the matrix symbol carries subscripts and 
superscripts, so we will use a special notation for Taylor/Laurent expansion coefficients of such matrix 
functions: if M|^(r) is such a matrix function and tq is a point about which this matrix is to be expanded, 
we write the Taylor expansion in the form 



OO 

En 



'M'jir-To)'. 



(1.5) 



We also use this notation in the case that tq — oo with obvious modifications. 

1.2. Assumptions and definition of fiuxon condensates. We study the Cauchy initial-value problem 



(1.1) under exactly the same assumptions used in our earlier work [7]: 
Assumption 1.1. The initial conditions for (1.1) satisfy F{x) = 0. 



Assumption 1 1 . 1| asserts that the initial data is of pure- impulse type. This is important because it implies 
that the direct scattering problem for the solution of the Cauchy problem by the inverse-scattering transform 
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Figure 1.3. Same as Figure [7T7| but with e = 0.046875. 




Figure 1.4. The initial data 0) — and ut{x,0) = G(a;) with G{0) < —2 plotted 
parametrically in the {u,UT)-plane (blue). Orbits of the simple pendulum are shown with 
gray curves. Orbits corresponding to nearby values of x (red curves) can diverge after a 
descent to a neighborhood of a saddle point in finite T, if x ^ ±a;crit — G^^{—2). 



reduces from the Faddeev-Takhtajan eigenvalue problem to the better-understood nonself adjoint Zakharov- 
Shabat eigenvalue problem. 

Assumption 1.2. The function G is a nonpositive function of Klaus- Shaw type, that is, G G L^(M)nC^(M) 
and G has a unique local (and global) minimum. 



As was shown by Klaus and Shaw in 121 , Assumption 1.2 provides a useful and important confinement 



property of the discrete spectrum of the nonself adjoint Zakharov-Shabat eigenvalue problem associated to 
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Figure 1.5. The region of the {x,t)-plane in which the fluxon condensate is analyzed. 



the potential G. This aUows us to formulate Riemann-Hilbert Problem |1.1| below on a system of contours 
very close to the unit circle and a transecting negative real interval, a set whose image under the function 
E{-) defined in (1.10) below is a segment of the imaginary axis. 

Assumption 1.3. The function G is even in x: G{—x) — G{x), placing the unique minimum of G at x — 0. 



Assumption 1.3 is admittedly less important (we believe that with appropriate modifications our results 
all hold true without it) but it allows for a substantial simplification of our analysis. We point out that 

has a well-defined inverse function G~^ : 



under Assumptions 
(G(0),0) ^ M+. 



1.2 



and 



1.3 



the function G restricted to 



Assumption 1.4. The function G is strictly increasing and real-analytic at each x > 0, and the positive 
and real-analytic function 



2G'(G-l(-^M))' 



< m < G{Oy 



(1.6) 



can be analytically continued to neighborhoods of m ~ and m — G{0)^ with '^{0) > and ^^(G(0)^) > 0. 

The analyticity of G for a; > and that of up to the endpoints of the interval [0, G(0)^] as guaranteed 
by Assumption 1.4 are both absolutely crucial to our method of analysis. It is the analyticity of G that 
implies that of '5 defined in (1.8) below, and hence of 6*0 defined by (1.10 1 and ( 1.14[ ), and these are used to 



convert a "primordial" Riemann-Hilbert problem of inverse scattering that involves a large number (inversely 
proportional to e > 0) of pole singularities into the simpler Riemann-Hilbert Problem |1.1| to be formulated 
below. The latter problem has no poles, but only jump discontinuities along contours, and hence is amenable 
to the Deift-Zhou steepest descent technique [16 of rigorous asymptotic analysis. 



Assumption 1.5. The small number e lies in the infinite sequence 

e = iV = l,2,3,. 



(1.7) 



Assumption 1.5 is important because it minimizes the effect of spectral singularities, events occurring 
infinitely often as e 4^ at which discrete eigenvalues are born from the continuous spectrum. When spectral 
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singularities occur, the reflection coefficient has poles in the continuous spectrum, and without Assump- 
tion 1.5 (or some suitable approximation thereof) the reflection coefficient cannot be neglected uniformly 



on the continuous spectrum. Assumption |1 .5| is important because our approach is based on neglecting the 



reflection coefficient entirely. (It can be shown to be small in the semiclassical limit without Assumption 1.5 
except near points where spectral singularities can occur.) 

Assumption 1.6. The function G{x) satisfies G{0) < —2. 

It is this last assumption that ensures that there exist exactly two points x — ixcrit at which the initial 
data curve {F,G) = (0, G'(x)) crosses the separatrix of the simple pendulum phase portrait. Therefore, 
Assumption |1.6| guarantees that the phenomenon we wish to study in this paper actually occurs for the 
initial data in question. 



Given initial data for the Cauchy problem (1.1 1 satisfying these assumptions, we construct a sequence of 
exact solutions u = un{x, t) of the sine-Gordon equation e^Ut* — e^M^a; +sin(u) = for e = ejv, N = 1,2,3,... 
called the fluxon condensate associated with the given initial data. While UN{x,t) is an exact solution of 
the partial differential equation for each iV, it generally does not satisfy exactly the given initial conditions. 
However, it has been proved [71 Corollary 1.1] that uj^{x,0) = ©(cat) holds modulo 47r and €NUN^t{x,0) = 
G{x) + 0{e]\[) where the error estimates are valid pointwise for x ^ and |x| ^ a;crit, and uniformly on 
compact subsets of the set of pointwise accuracy. We strengthen this convergence to include the points x 
near where |a;| = a;crit in Theorem |1.1| below. 

The fluxon condensate is constructed as follows. First one defines a function ^'(A) by setting 

^iX):^ir ' ^Gisr-v-^ds, A=^, 0<i;<-G(0). (1.8) 

Note that 5* is a strictly decreasing function of v where defined. This fact allows us to define a sequence of 
numbers {X^n k}k=o solving the equation 

vI/(A^,fc) = Trejv (^fc + ^) , A; = 0,l,2,...,iV-l. (1.9) 
(These positive imaginary numbers are approximate eigenvalues for the Zakharov-Shabat eigenvalue problem 



associated with the Klaus-Shaw potential G, and (1.9) is a kind of Bohr-Sommerfeld quantization rule for 
that nonself adjoint problem.) Setting 

i 



E[w) \ 



{-wY'^ + {-w)-^/"^]^ and £>(«;):==- [(-u;)i/2 - (-u;)-i/2j , | arg(-'u;)| < tt, (1.10) 
we define 

Q{w]x,t) -.^ E{w)x + D{w)t, I arg(-w)| < TT, (1.11) 

and 

■■= n J |arg(-HI < (1.12) 

This function is meromorphic and has no zeros where defined. It has 2A^ poles (counted with multiplicity), 
the set of which we denote as Pn- Generically (with respect to deformations of G) the poles are all simple, 
and there are 2Ns of them in complex-conjugate pairs on the unit circle in the w-plane, and 2A'^k of them 
in pairs on the negative real axis in involution with respect to the map w 1/w. AsA^fcxa (or ^ 0) the 
poles accumulate on the whole unit circle and the interval [a, b] , where 

a:--^(VG(0)2-4-G(0))', b := ( ^6(0)^ - 4 + G(0)) ' ^ ^. (1.13) 

Clearly, both a and b are independent of N, and o<— l<b<0. Were Assumption |1.6| not satisfied, 
we would have A^k — 0, and the poles of IIjv would accumulate only in a complex-conjugate symmetric 
incomplete arc of the unit circle containing the point w = 1. In nongeneric cases it can happen that there is 
a double pole aX w = —1. We consider the accumulation locus of poles. Poo := [a, b] U 5"^ to be an oriented 
contour, with orientation of the two intervals (a, —1) and (—1, b) both toward w = —1, and with orientation 
of the upper and lower semicircles from w = —1 toward w = 1. See Figure [L6j We also use the notation 

eo{w):=^{E{w)). (1.14) 

7 



Then, define a function by a Cauchy integral: 



L{w) 



and then set 



Finally, writing 



YNiw) 



^ r eo{y) dy 



u.eC\(PooUM+), 



1 



m :=-(L+(0 + i-(0), 



w e c \ (Poo u : 



(1.15) 

(1.16) 
(1.17) 



a function that turns out to have a well-defined analytic continuation to a full neighborhood of the self- 
intersection point of w = — 1 of the contour Poo (upon continuation from any of the four intersecting arcs), 
we set _ 

Tjv(?) :=2nAr(e)cos(e^i0o(e))e-^^«^/^", C e Poo- (1.18) 
It can be shown W, Proposition 3.1] that Ytq(w) = 1 -|- C(e7v) holds uniformly on compact subsets of the 
open domain of definition, and that YAr(w) = 0(1) holds if kj a or w — b nontangentially to the real 
axis. Similarly, Tn{^) extends from Poo to an analytic function Tn{w) on a neighborhood of it; ~ —1, and 
Tn{w) = 1 -|- 0{£n) holds uniformly on this neighborhood as well as on Poo, as long as w is bounded away 
from w = a and w = b; however the estimate Tn{^) =0(1) holds uniformly for ^ € Poo- 



Consider the contour En illustrated in Figure 1.6 Let giw) denote a function analytic in the domain 



R+\/ y 






Ri\I 





■ OO 



Figure 1.6. The contour En of discontinuity for the sectionally analytic function N(u'). 
With the exception of the two components of M+ \ / which are oriented left-to-right, the 
contour arcs are oriented with the regions labeled "-\- " on the left, and with the regions labeled 



"—" on the right. The image o/Sn\K+ under w i— >■ E{w) with E defined by (1.10) consists 
of the four straight line segments: (i) "^{E} = -5,Q < '^{E} < -G(0)/4, (ii) 5R{P} = 0, 
< "^{E} < -G(0)/4, (iii) 5R{P} = (5, < Q{E} < -G(0)/4, and (iv) |3?{S}| < 5, 
= — G(0)/4. Here 6 > is a sufficiently small number. The self-intersection point 
w — —1 is the critical point of E. 

w E C \ (Poo U M+), satisfying the symmetry g{w*) = g{w)* and the conditions that ^(0) = and g{w) — > 
as w ^ OO. We assume also that g takes well-defined boundary values on Poo U M+ in the classical sense 
(Holder continuity up to the boundary), and that the boundary values satisfy .g+(^) -f .g-(C) = for ^ € . 
Set 

0(0 ■■= -^{g+iO - 9-iO) and 0(e;x,t) :=2*Q(C;a;,i) + i(0-.9+(0-5-(e), C ^ Poo- (1-19) 
Then, to determine the fiuxon condensate, one solves the following Riemann-Hilbert problem. 

Riemann-Hilbert Problem 1.1 (Basic Problem of Inverse Scattering). Given g as above, find a 2 x 2 
matrix function N(7«) = N^(it;; x, t) of the complex variable w with the following properties: 



Analyticity: N(zi;) is analytic for w G C \ En- 

Jump condition: Nlw) is Holder continuous up to En- The boundary values N±(^) taken by 
N(w) on the various arcs o/En are required to satisfy the following jump conditions 

N+(o -tT2N_(e)a2, 

1 + (,^{Sa+(i)-eo-{i))/<^N iy^_(^)e-i«o-(«)/eiVe-F(C:2;,t)/ejv" 



N+(0 = a2N_(0a2 
where x, t) 



■.^2iQ+{tx,t) 



N+(0 = N_(0 



jeo+(?)/eNgF(C:x,t)/e„ 

g-je(C)/£iv 



1 



(1.20) 

Ce/, (1.21) 





e»e(?)/«N 



and, 



N+(e)=N_(0 



_jy^(^)e(2«Q(?;2:,t)+i(«)±ieo(C)-23(c))/ejv 
Normalization: The following normalization condition holds: 

lim N(w) I, 

liJ— >-oo 

where the limit is uniform with respect to angle for \ arg(— < tt. 



? e Ej 



(1.22) 

(1.23) 
(1.24) 



It turns out that for each choice of the function g, each N — 1,2,3,..., and each (x,t) e there 
exists a unique solution of this Riemann-Hilbert problem. Moreover, the product M(w) = 'M.n{w; x,t) := 
N^(?«; X, i)e^*-"'-''^^/^" does not depend on the choice of the function g. The solution N^(w; x, t) has conver- 
gent expansions near w = and w = cxi of the form 



and 



N%{w;x,t) =1- 



N%{w;x,t) 



^'Nf^]{x,t){—w) for jwl sufficiently large. 



k=0 



lN%]{x,t){-wf/^, for \w\ sufficiently smah. 



(1.25) 



(1.26) 



Definition 1.1 (Fluxon condensates). Given a function G satisfying A ssumptions \l . i] through 1.6 the fluxon 
condensate associated with G is the family of functions {uN{x,t)}^^-^ defined modulo Att by the equations 



cos.{\uN{x,t)) -.^ [oN^](a;,t) 



and 



11 



sin(iuAr(x, t)) 



21 



(1.27) 



(1.28) 



Note that these are independent of g because g{Q) = 0. 

While it is possible to extract formulae for derivatives of UN{x,t) with respect to x and t by the chain 
rule, it is preferable to have formulae that do not require differentiation of N^(i(;; a;, <:), as this will require 
control of derivatives of error terms. However, it can be shown also that the following formula holds: 

^""^ (x, t) = [HN^^] (:r, t)] + [[X^] (x, O-M^N^] (x, , 

and this does not require differentiation of N^(w; x, t) with respect to t (it also turns out to be independent 
of choice of g). Each function UN{x,t) of the fluxon condensate is an exact solution of the sine-Gordon 
equation with e — e^r, and ujv(— a;,i) — ui\[{x,t). 

While Riemann-Hilbert Problem |1.1| is the most convenient starting point for our local analysis near 
criticality, it is not the most fundamental representation of UN{x,t). Indeed, UN{x,t) is a reflectionless 
potential for the direct scattering problem of the Lax pair for the sine-Gordon equation, and this means that 
it can be obtained from a purely "discrete" Riemann-Hilbert problem whose solution is a rational function 
on the Riemann surface of y"^ — —w, with poles on both sheets over the points of P/y- Since the number 
of poles is increasing with N, some preparations are required to recast the problem in a suitable form for 
addressing the limit N ^ oo. These preparations are detailed in [7^, and they take the form of a sequence 
of explicit transformations resulting in the equivalent Riemann-Hilbert problem |1.1[ In general, a number 
of choices are made along the way because the transformations that are required to enable the subsequent 
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asymptotic analysis depend on {x,t); however for x « XcHt > and t small only the simplest of the choices 
detailed in [7] are required. For readers familiar with the terminology of that paper, we are assuming that 
A = 0, and correspondingly, YAr(w) is the function called ¥^^{10) in [7j while Tn{w) is the function called 
T"^ (w) in YJ. Also, for our purposes we need make no distinction between the contour called E^UE'^ (which 
would just be E^ in the case that A = 0) and the contour Poo, and this implies that the set Z C C that lies 
between these two contours is empty. Finally, to derive Riemann-Hilbert Problem from the results of [7] 
we used the facts that for ^ G /, Q+{^;x,t) + Q^{(;x,t) L+(0 + -^-(C) = 0. 

The solution of the sine-Gordon equation by inverse-scattering methods is a subject with a long history. 
When the sine-Gordon equation is written in terms of characteristic coordinates a; ± t it fits naturally into 
the hierarchy of the Ablowitz-Kaup-Newell-Segur or Zakharov-Shabat scattering problem. The characteristic 
Cauchy problem for the sine-Gordon equation was integrated in this framework by Ablowitz, Kaup, Newell, 
and Segur [I] and by Zakharov, Takhtajan, and Faddeev [21]. The more physically-relevant problem of the 



Cauchy problem in laboratory coordinates as posed in (1.1) required new methodology, and the solution 
of this Cauchy problem by the inverse scattering method was first outlined by Kaup |2Q]. An account 
of the solution of the Cauchy problem in laboratory coordinates is given in the textbook of Faddeev and 
Takhtajan |17j . Some further analytical details needed to make the theory completely rigorous were supplied 
by Zhou and later by Cheng, Venakides, and Zhou [5] . A self-contained account of the Riemann-Hilbert 



formulation of the inverse-scattering solution of the Cauchy problem ( 1.1 ) assuming only that at each instant 
of time the solution has L^-Sobolev regularity can be found in our paper [Q, Appendix A], and a direct 
proof that the sine-Gordon equation preserves this degree of regularity if it is initially present is given in [SJ 
Appendix B] . In our recent paper [J| we found that to describe the semiclassical limit for the Cauchy problem 



(1.1) it is useful to reformulate the Riemann-Hilbert problem in the complex plane of a square root of the 
spectral variable used in [SI Appendix A]; this makes it easier to express the asymptotic solutions in terms 
of Riemann theta functions of the lowest possible genus. As we view the current paper as a continuation of 
our work in we use the same formulation here. 

1.3. Statement of results. Let Xcrit be defined by 

^Ecrit :=G"'(-2) >0, (1.29) 

and define a positive constant v by 

> 0. (1.30) 



12G'(Xcrit) 



Also, set 

Aa: := a; - a;crit- (1-31) 
All of our results concern the asymptotic behavior of the fluxon condensate un {x, t) in the small region near 
criticality where t = 0{€n log(e^^)) and Ax = ©(e^^) as shown in Figure 



1.5 



Our first result is concerned with the relevance of the fluxon condensate associated with G{x) to the 
Cauchy initial-value problem. 

Theorem 1.1 (Initial accuracy of fluxon condensates). Suppose that x ± a^crit = ^(c^'^)- Then uniformly 
for such x, 

UN{x,Q)=0{e^J;^) (mod47r) and tN-^{x,Q) ^ G{x) + 0{e^J;^). (1.32) 



This result extends that of [71 Corollary 1.1] to suitable neighborhoods of the points x — ±Xc 



Theorem 1.2 (Main approximation theo rem). There exist multiscale approximating functions C(a;,i;ejv) 
and S{x,t]eN) (defined in detail in t'7.3, and depending on initial data only through the constant v > 0) 
such that 



coa{luN{x, t)) = C{x, t; eat) + 0{e]^^) 
sin(iuAr(x, t)) = S{x, t; eAr) -|- 0{e]l^) 
with the error terms being uniform for Ax ~ 0{e^J^) and t ~ 0{eN log(e^^)). 



(1.33) 
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In fact, we will really show that the error term is significantly smaller, namely 0{e^ ), over most of the 
small region of the (x, i)-plane where the above Theorem provides an asymptotic description of the dispersive 
breakup of the pendulum separatrix. Now, the multiscale model provided by Theorem 1 1 . 2| serves to establish 
universality of the behavior near the critical point, but the formulae for C{x, t; e^r) and S{x, t; e^) are rather 
complicated, so it is useful to give some more detailed information by focusing on smaller parts of the 
{x, t)-plane near criticality. 

To render our results in a more elementary fashion, we need to first define a certain hierarchy of rational 
functions. First define 

llo(z) 1 and Vo(z) --z. (1.34) 

6 

Then define {(Urm^m)}mei. by the forward recursion 

U^+,iz):^~lzU^iz)~^^f^ + lu';,{z) and V„+i(z) (1.35) 
and the backward recursion 

U^-i{z):=-^ and V™_i(z) ^V:;(z) - ^^f^ - ^zV™(z). (1.36) 

Vm[Z) 2 Vm{z) 6 

Up to constant factors, Um{z) and Vm{z) are ratios of consecutive Yablonskii-Voroh'ev polynomials (see 
[13l[T4]) and logarithmic derivatives oiUm{z) and Vm{z) are the unique |22] rational solutions of the inho- 
mogeneous Painleve-II equations. These observations are not necessary for us to state our results, and we will 
make further comments later at an appropriate point in the paper. The following proposition characterizes 
the behavior of these rational functions for large z. 



Proposition 1.1. For each m £ Z, 

l — m 

as z — > oo. In particular, sgn(Um{z)) — sgn(Vm(z)) = 1 for sufficiently large negative z while sgn(Z^,„(z)) = 



(l + 0(z-i)) and V„(z)=(--) {l + 0{z-^)) (1.37) 

rticular, sgn(Um{z)) = sgn(Vm(z)) = 
— sgn(V„i(z)) — (—1)™ for sufficiently large positive z. 

Proof. It is obvious that the recursions ( 1.35 1 and ( 1.36 ) preserve rationality of the input functions (Um, V„i), 
so for each to G Z, Lim is a rational function of z. Therefore, Um has a Laurent expansion for large z € C of 
the form 

U^iz) = c^z''^il + 0{z-^)) (1.38) 
where Cm is a complex constant and km G Z is an exponent to be determined, and moreover this expansion 
is differentiable any number of times with respect to z. In particular, it follows that 

^^f^ = klcmz'"^-\l + 0{z-')) and Ull^iz) ^ km{k,n ~ l)c,nz'''--\l + Oiz'')) (1.39) 

as z — ^ oo. This shows that the final two terms on the right-hand side of the formula for Um+i{z) given in 
(1.35) are subdominant compared to the term —zUmiz)/6, and therefore, 

Um+liz)^~Umz'''-il+0{z-'))^-^Umiz){l+Oiz-')), Z-^OO. (1.40) 

o 6 
This formula gives the recurrence relation for the large-z asymptotics of Um- The asymptotic formula for 
Um{z) given in (1.37) follows by solving the recurrence with the base case Uq{z) — 1, and that for Vm(z) 
follows from the identity Vm{z) — l/Um-i(z). □ 

Theorem 1.3 (Superluminal kink asympt 
\t- |TOeArlog(e^^)| < |eArlog(e^^). Then 



Theorem 1.3 (Superluminal kink asymptotics). Fix an integer m, and suppose that Ax — C'(e^'^) while 



where the phase is 



cos{^UNix,t)) = {-ly' sgn{Um{z))sech{TK) + Ecosix,t;eN) 
sin(iujv(a;,i)) = (-ir+'tanh(TK) + £;sin(x,t;e7v), 

t Ai/i/A 
Tk — -2TOlog I I +\og\Umiz)l (1.42) 



and a rescaled spatial coordinate is given by 



The error terms satisfy 



lim Ecosix,t;eN) = 1™ Esinix,t;eN) = 



(1.43) 



(1.44) 



unless (i) both e^^^/^e"*/^" = ^(4^^) and \z - 3f{Um)\ = 0{e]^^) or (ii) both e^™/^e*/<^" = ^(4^^) ^'^'d 
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\z — i^ilArn)] — C'(ejv )■ Moreover, given any interval on which \og\Um{z)\ is a bounded function, 

we have 

i;cos(x,t;ew) =0(4') and E,,^{x,t-eM)^0{e^jl'') (1.45) 
holding uniformly for z G and \t — |mejv log(e]^"'^)| < |eAr log(e^"'^). 

Wherever the error terms vanish in the limit, it follows that cos{ui^{x,t)) — > 2sech^(rK) — 1 and that 

2tt by 

(1.46) 



sai{uis[{x,t)) — > — 2(7 sech(TK) tanh(TK), where a := sgn(t/,„(z)). If we define a function u{T) modulo 27r by 

cos(u(T)) := 2sech^(r) - 1 

sin(u(r)) := -2(Tsech(r) tanh(r), 
then it is easy to check that u{T) is an X-independent solution of the unsealed sine-Gordon equation 



+ sin(7i) = 0. 



(1.47) 



This exact solution represents a superluminal (infinite velocity) kink with unit magnitude topological charge 




T 



-0.5 







0.5 



Figure 1.7. (Superluminal antikink). The cosine (left) and sine (right) ofu(X,T) = u(T) 
defined by (1.46) with a = I. For a superluminal kink (a = — Ij, the plots are simply 



reflected through the horizontal X-axis. 



cr; if (T = 1 (respectively a = —1) then u decreases (respectively increases) by 27r as T varies from T 



to r = +OD. Sometimes one refers to u(T) as a kink in the case cr = —1 and as an antikink in the case a — 1. 
Another way to describe u{T) is to say that it is a solution of the simple pendulum equation u" + sin(it) = 
that is homoclinic to the unstable equilibrium of a stationary inverted pendulum. See Figure [TTtI 



The dominant part of the phase variable Tk defined by (1.421 is indeed a recentering and rescaling by 
of t; however Tk also contains weak z-dependence through the function log |ZYm(2:)| where z is proportional 
to X — Xcrit- Therefore, the approximation of UN{x,t) described in Theorem 1.3 is not an exact kink, but 



rather is one that is slowly modulated in the direction parallel to the wavefront. Indeed, the center of the 
approximating kink (where the pendulum angle is zero) corresponds to Tk = 0, which is a curve in the 
(x, i)-plane that is a scaled and translated version of the graph of the function — log \Um{z)\- See Figure 1.8 
An additional important observation is that the period of each approximating kink is proportional to e^r, but 
the strip in which it lives (indexed by the integer to) corresponds to a time interval of length proportional to 
e^v log(e^ ); hence the kinks are in reality widely separated from each other, and therefore in 
domain lS.x = C'(e^'^) and t — 0{eN log(e^^)) covered by Theorem 
formula UN{x,t) w tt (mod 27r) holds, which in the context of the coupled pendulum interpretation of 
sine-Gordon means that "most" of the pendula are approximately in the unstable inverted configuration. 



most" of the 
it is fair to say that the approximate 
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Figure 1.8. The horizontal strips in the {z,t)-plane in each of which Theorem 1.3 provides 
asymptotics for U]\[{x,t), shown with alternating gray and white shading, with the center 
Tk = o/ each corresponding superluminal kink shown with a superimposed curve. Left: 
tjsi = 10^'"'. Right: — lO^''. In both cases Av^^^ — 1. For different values of v > 0, the 
curves can "march out" of the confining strips as \t\ increases, unless cat is small enough. 



We may now formulate some observations about the kink centered at Tk = that follow from Propo- 
sition [lT] Indeed, for each m G Z the following is true: for z sufficiently negative, a := sgn(Z^m(z)) = 1, 
while for z sufficiently positive, a — (—1)"*. This then implies that as z — >■ — oo we have a series of time- 
like antikinks (all of the same topological charge) consistent with nearly synchronous rotational motion of 
pendula. Similarly, as z — !■ +oo we have alternation between kink and antikink with no net topological 
charge consistent with nearly synchronous librational motion of pendula. These facts demonstrate that the 
asymptotic behavior near criticality matches appropriately with the established asymptotic formulae ij^ valid 
away from criticality that model uiq{x,t) as a modulated train of superluminal rotational traveling waves 
to the left of Xcnt and as a modulated train of superluminal librational traveling waves to the right of Xcnt . 
It also follows from Proposition |1.1| that if to > then Um{z) blows up as z — >■ ±oo while if to < then 
l^m{z) — > as z -> ±oo. This implies that for to > the corresponding kinks follow a logarithmic "frown" 
for large |z|, while for to < they instead follow a logarithmic "smile" for large |z|. These features can be 
seen in Figure |1.8| 

When the horizontal strips of the {x, t)-plane corresponding to different integral values of to are put 
together, one sees that the kink-like asymptotics given by Theorem |1.3| are valid throughout the region 
where Aa; = ©(e^/'^) and t = O{eN^0g{e^)) with the sole exception of small sub-regions near the top and 
bottom edges of each strip where the curve Tk = tries to exit the strip. When the strips are put together 
there are obvious mismatches of the curves from neighboring strips in these small sub-regions as can be seen 
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from the left-hand plot in Figure 1.8 Our final result corrects the kink asymptotics in these regions and 
therefore removes the mismatches. 



Theorem 1.4 (Grazing kink coUisional asymptotics). Fix an integer m and let zq denote any of the 
(necessarily simple) real zeros of U„i-i{z) (equivalently a simple zero ofVm+i{z)). Given a sufficiently 



small positive number /i > 0, suppose that \t — (f^^ ^ |)£-/v log(ejv^)| < ^Eat log(e^"'^), and that \z 
fj,e]^^ exp{\t — (|m— ^)eN log(e^^)\/ (2ej^)) , where z is defined by (1.431. Then 



< 



2XGsech(rG) 
1 + Xlscch^{Tc] 



+ Ecosix,t;eN) 



l-X|sech'(TG) 
1 + X^sech2(TG) 



(1.48) 



where 



Xc 



(-N 



1/3 



Ax 



2^.1/3^2/3 



and Tc 



t 
CAT 



(2m - 1) log 



4^,1/3 



1/3 



The error terms satisfy 



lim Ec, 



s{x,t;eN) = lim Es- 



^{x,t;eN) = 



logK_i(zo)|. (1.49) 



(1.50) 



whenever \z — Zq\ <^ I as ejv i (which excludes only the four extreme corners of the hourglass-shaped region 
under consideration). Moreover, we have 



Ecos{x,t;eN) 



holding uniformly for \z — Zq | 



o(ey^^ 



and \t 



and Esinix, t;eN) 



(1.51) 



(frn- i)ejvlog(e^^)| < iejvlog(e/^ 



Wherever the error terms vanish in the limit, it follows from (1.48) that cos{uN{x,t)) — > cos(u(Xg, Tg)) 
and sm{uN{x,t)) — > sin(M(XG, Ig)), where with k := sgn(W^j_]^(zo)), a function u{X,T) is defined modulo 
2tt by 



cos(m(X,T)) := 



sin(u(X,T)) := 



8X2sech^(T) 



(1 

4KXsech(r) 



- 1 



X2sech2(T))2 

l-X2sech^(T) 



(1.52) 



(1 + X2sech2(r))2' 

Then it is again easy to confirm that u{X,T) is an exact solution of the unsealed sine-Gordon equation 
(1.47). Plots of cos(u(X, T)) and sin(u(X, T)) are shown in Figure 



1.9 



This is a particular solution of the 
sine-Gordon equation that corresponds to boundary conditions of the form m — >■ tt (mod 2tt) as T — ^ ±oo. 
In the proper version of scattering theory corresponding to these boundary conditions, the solution at hand 
is a reflectionless potential associated to a double eigenvalue. In the context of the Riemann-Hilbert problem 
of inverse scattering such an object is encoded as a double pole of the matrix unknown. Such solutions were 
noted in the earliest days of inverse scattering by Zakharov and Shabat [53] in their study of the focusing 
nonlinear Schrodinger (NLS) equation. They pointed out that such solutions describe grazing collisions of 
solitons. The trajectories of the solitons emerging from the interaction region are not asymptotically straight 
lines, but rather are logarithmic, with relative velocities of the solitons tending to zero. In the context of 
the sine-Gordon equation we have instead a grazing collision of superluminal kinks. These structures serve 
to smooth out the mismatches of curves shown in Figure [L8] at horizontal strip boundaries. In the language 
of matched asymptotic expansions, they function as internal transition layers of "corner" type. 



The regions of validity of the asymptotic formulae given in Theorems |1.3| and 1.4 are compared in Fig- 
ure |1.10[ which illustrates the complementary nature of these two results. It is clear from this figure that 
Theorems |1.3| and |1.4| both provide simple asymptotic formulae for the same fields in overlapping regions. 
In the overlap domains, formulae of kink type provided by Theorem |1.3[ and of grazing kink collision type 
provided by Theorem 1 1 .4| are simultaneously valid. This is actually a consequence of Theorem |1.2[ but it may 
also be checked directly, by writing both types of formulae in terms of common spatiotemporal independent 
variables. 
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Figure 1.9. (Grazing collision of superluminal kinks). The cosine (left) and sine (right) 
of u{X,T) defined by (1.52) with k = 1. For k = —1 the plots are simply reflected through 
the vertical T-axis. 



One way to characterize the asymptotic formulae giving the universal form of the dispersive breakup of the 
simple pendulum separatrix under the sine-Gordon equation is to say that the wave pattern consists of waves 
of elementary "solitonic" forms that are spatiotemporally arranged according to solutions of certain nonlinear 
differential equations of Painlcve type. To our knowledge, this type of phenomenon was first observed quite 
recently in a paper by Bertola and Tovbis ^ in which all of the local maxima of the modulus of the solution 
to the focusing NLS equation in the semiclassical limit near the onset of oscillatory behavior (the point of so- 
called elliptic umbilic catastrophe of the approximating elliptic quasilinear Whitham modulational system) 
are individually modeled by the same exact solution of the focusing NLS equation (the rational or Peregrine 
breather solution) but the locations of the maxima are far more transcendental, being determined by the 
location of poles in the complex plane of the so-called tritronquee solution of the Painleve-I equation. In 
the current context, the role of the tritronquee solution is played by the rational functions {Um}mi£Z solving 
the Painleve-II system, and interestingly, we require an infinite number of different Painlcve functions to 
describe the full wave pattern. 



2. Choice of g{w) 

The use of a so-called "g-function" to precondition a matrix Riemann-Hilbert problem for subsequent 
asymptotic analysis by the Deift-Zhou steepest-descent method is by now a standard tool, having first been 
introduced by Deift, Venakides, and Zhou [15] in their analysis of the Korteweg-de Vries (KdV) equation in 
the small-dispersion limit. Recall that the scalar function g appearing in the formulation of Riemann-Hilbert 
Problem |1.1 1 is subject to the basic requirements of 

• Analyticity: g is analytic in the domain C \ (Poo U 

• Boundary behavior: g takes well-defined boundary values on Poo U R+ in the classical sense, and the 
sum of the boundary values vanishes on and, 

• Normalization: g{0) = g{oo) = 0. 

We can use any such function to formulate the basic Riemann-Hilbert problem for N(i(;) — N^(?i;; x, t), 
but it is well-known that to study the limit — )■ oo it is important that g be chosen appropriately. In 
[7] it is explained how g should be chosen given {x,t) e M.'^ to calculate the limit N ^ oo with {x,t) 
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Figure 1.10. Theorem] 1.3\provides asymptotics in horizontal strips in the {z,t)-plane with 

1 /3 

the exception of the small rectangles shown in green, and with optimal accuracy of 0{€j^ ) 
in the suhregions shown in blue. Theorem I.4 provides asymptotics in hourglass-shaped 
regions with the exception of where they overlap the blue regions ( where the error estimate 

1 /3 

of Theorem 1.3 is optimally small) and with optimal accuracy of 0(ej^ ) in the included 



green regions (where the error estimate of Theorem 1.3 fails to be small) and brown regions. 



fixed. The chosen function g{w) = g{w;x,t) has the property that its boundary values satisfy a kind of 
equiUbrium condition (an analogue of (2.2 1 below) on a contour (called /3 in [7]) whose topology changes at 
criticality: in the simplest (genus one) case, near {x,t) — (xcrit,0) the equilibrium contour is either a nearly 
circular arc through w ^ I with complex-conjugate endpoints near w = — 1 (this corresponds to modulated 
superluminal librational waves) or the union of a nearly circular closed contour with a transecting small 
interval of the real axis with both endpoints near w ~ —I (this corresponds to modulated superluminal 
rotational waves). Exactly at criticality the two endpoints coalesce at w = — 1 and the equilibrium contour 
becomes the unit circle. It turns out that proximity of endpoints of this contour is an obstruction to proving 
uniform asymptotics for (x, t) near, but not exactly at, criticality. 

Here, since we want to allow {x,t) to approach the critical point (XditjO) at some rate depending on 
iV, we need to use a different function than g(w) = g{w]X,t) defined in |7]. We will replace g{w) with 
q{w) — Q{w;x,t) whose construction and properties are described below. It is important that while for 
general x and t, the functions g{w) and q{w) are different, at criticality they coincide (they satisfy the same 
conditions, which can be shown to uniquely determine the solution). 

The function q{w) is required to satisfy the following conditions: 

• q{w) is analytic if \w\ 7^ 1 and w ^ M+, and it takes continuous and bounded boundary values on 
the unit circle = 1 and on the positive half-line IR-|_. 

• q(w) obeys the following jump conditions: 

0+(C) + 0-(C)-O, ^eM+. (2.1) 

2zQ(C;a;,t) + L(e)-0+(O-fl-(e) = O, 1^1 = 1. (2.2) 

• q{w) has the following values: 

0(0) = 0. (2.3) 

lim q{w) = 0. (2.4) 



16 



To obtain a formula for 5{w), we first write q in the form 

g{w)^h{w) + t){w), (2.5) 

where the function L is defined by (1.151. Therefore the above conditions on g are translated into equivalent 
conditions on (): 

• i){'w) is analytic if \w\ ^ 1, w ^ M+, and w ^ [a, b], and it takes continuous boundary values on the 
unit circle \w\ — 1, on the positive half- line M+, and on the interval [o, b]. 

• ()('!«) obeys the following jump conditions: 

fl+(O + f)-(O = 0, ^€M+. (2.6) 

2iQ{C,x,t)-l)+{O-^-iO^0: iei = l- (2.7) 

l)+{0-^-{0 = -UL+{0-L-{0) = -^eo{0, we[a,b]. (2.8) 
In the latter condition, the two complementary segments (a, —1) and (—1, b) are taken to be oppo- 
sitely oriented toward ^ = — 1. 

• i){'w) has the following values: 

f)(0)=0. (2.9) 

lim l){w) = 0. (2.10) 

Finally, we write Ijiw) in terms of a new equivalent unknown m{w) by the piecewise substitution: 

^^^) ■= \ I ^-l/2 ( \ 11^1 2.11) 
!-(-!«) ^/"^m(w), \w\ < 1. 

The conditions satisfied by m{w) are then the following: 

• m{w) is analytic if \w\ ^ 1 and w ^ [a, b], and it takes continuous boundary values on the unit circle 
\w\ = 1 and the interval [a, b]. 

• m{w) obeys the following jump conditions: 

m+(0-m_(0 = -2^(-0'/'Q(C;a:,^), 1^1 = 1, (2-12) 
where the orientation of the circle is in the positive (counterclockwise) sense, and 

m+(0-m-(0 = -*(-0'/'^o(e), ee[a,b], (2.13) 

where now the interval [a, b] is taken to be oriented left-to- right. Note that (— ^)^/^ > for a < ^ < 
b < 0, and that the right-hand side of (2.12) is a single- valued function on the circle: 

-2z(-6'/'Q(^;^,t)-|(i-0-^(i + 0, iei = i. (2.14) 

• m(w) is uniformly bounded, and has the following value: 

m(0) = 0. (2.15) 
We may now express m{w) in terms of a Cauchy integral using the Plemelj formula: 

, , ^ 1 / ~2»(-g)V^Q(e;x,t) , 1 ^thO'/'9o{0 ,„ 

where the constant mo is the most general entire function we can add to the Cauchy integrals consistent with 
the uniform boundedness of m(w). The integral over the positively-oriented unit circle may be evaluated in 
closed form: 

m(^)^mo + ^^ e 1 '^-^ ^'-''^ 

0, kl>l. 



We choose the value of the constant mo to enforce the condition (2.15): 
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It follows that 



m{w) = 



X + t 



-w, 



\w\ < 1 



\w\ > 1. 



(2.19) 



By writing g in terms of I) (by (2.5)) and then writing [} in terms of m (by (2.11 1), this completes the 
construction of Qiw). 



3. Steepest Descent, the Outer Model Problem, and its Solution 



According to (2.2 1, with N(ti;) = N^(w;x,t) and with q{w) — g{w;x,t) defined as in ^ a condition of 
"equilibrium type" holds on the unit circle \w\ = 1, regardless of the values of x and t. The next step is to 
"open a lens" about the whole circle. The lens consists of four disjoint open sets separated by the unit circle 
and the interval / C 



as shown with shading in Figure 3.1 



The union of the two components lying on 
the left (respectively right) of the unit circle with its orientation will be denoted A+ (respectively A~). We 
now define a new unknown 0{w) equivalent to N(w) by the substitution 



N(w)rAr(w)-''^/2 



0(w) 



1 =Fie 




-[2iQ(w;x,t)+L{w}Tieo{w)-2g{w)]/eN 



€ A± 



w e C \ (En U A+ U A- 



(3-1) 

Here, the square root Tn{'wY^'^ is well-defined as the principal branch, since the uniform approximation 
Tn{w) ~ 1 holds when N is large. (The definition (3.1 ) coincides with that of 0{w) in terms of N(u;) given 
in [7] in the case that A = and with g replaced by g.) The jump contour for 0{w) is the same as that for 
N(w) but augmented by four arcs emanating from w = —1 representing the "outer" boundaries of the lens 
halves A* as shown in Figure 3.1 




Figure 3.1. The contour of discontinuity for the sectionally analytic function 0{w), with 
the lens halves A^ indicated with shading. Arcs of the jump contour also present in the 
jump contour for TSI (w) (see Figure 



1.6) retain that orientation, and the lens boundaries are 



oriented as indicated. 



It is a consequence of the equilibrium condition (2.2) satisfied by the boundary values taken by q{w) on 
the unit circle that 0{w) satisfies the piecewise-constant jump conditions 



O+(0 = O_(0(T*^i), ICI = i, ±GU}>o. 
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(3.2) 



Since the equilibrium contour is a closed curve, we may remove these discontinuities from the problem by 
making another explicit transformation: 



P{w) 



'OH, 

0{w){-i<Ji) 
0{w)iai, 



\w\ > 1 

\w\ < 1 and > 

< 1 and < 0. 



(3.3) 



Unlike 0{w), the matrix function P(w) extends continuously to the unit circle, but the jump conditions it 
satisfies within the unit disk where it differs from 0{w) are altered somewhat from those of 0{w), including 
a new jump of the simple form P+(^) = — P-(^) on the segment b < ^ < 0. The contour of discontinuity 



for P{w) is illustrated in Figure 3.2 




Figure 3.2. The contour of discontinuity for the sectionally analytic function P(u'). Con- 
tour arcs also contained in the jump contour for 0{w) (see Figures 1.6 and 3.1) retain that 
orientation, and otherwise (that is, in the interval (b,0)j the orientation is as indicated. 



The jump discontinuities of P{w) will turn out to be negligible except in a small neighborhood of w — —I 
and along the ray w > —1, and P{w) is a matrix tending to the identity as w — >■ oo. These considerations 
lead us to propose a model Riemann-Hilbert problem whose solution we expect to approximate P{w) away 
from w = —1: 

Riemann-Hilbert Problem 3.1 (Outer model problem near criticality). Let w^, w — 1 6e a real parameter. 
Find a 2 X 2 matrix P°"*(w) with the following properties. 

Analyticity: P°"*(w) is an analytic function of w for w G C \ [w^,+oo). Holder continuous up 
to the jump interval [w^,+oo) with the exception of an arbitrarily small neighborhood of the point 
w — w^. In a neighborhood of w ^ w^,, the elements o/P°^'(w) are bounded by an unspecified power 
o/ |w - 

Jump condition: The boundary values taken byP°^^{w) along {w^,0) andR^ satisfy the following 
jump conditions: 

P°f{0 = -P°-HO, w,<^<0, (3.4) 

and 

Pr(0 = ^2P™*(0f^2, ^eR+. (3.5) 

As these jump conditions are both involutive, orientation of the jump contours is irrelevant. 
Normalization: The following normalization conditions hold: 

limP°"*(w)=I and det(P°"'(w)) ee 1. (3.6) 
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There are many solutions of this Riemann-Hilbert problem, a fact that can be traced partly to the 
unspecified power-law rate of growth admitted as w — >■ . The complete variety of solutions is not directly 
relevant for us; we will simply select a family of particular solutions in the class of diagonal matrices. Indeed, 
for each m G Z, we have the solution 



a/2 



(l-2m)<T3 



W 



(m-l/2)o-3 



(3.7) 



Note that if m is held fixed, then P^*(w) and its inverse are bounded when w is bounded away from w^. 
Later (see (4.52|), we will let depend weakly on x and t in such a way that — )• —1 as t — )• and 
X — >■ Xcrit- Clearly, any bound for P°"*(w) or its inverse valid for w bounded away from —1 will hold 
uniformly with respect to x and t near criticality, and of course P^'(?i') is independent of e^r. 



4. Inner Model Problem Valid near w = -1 

4.1. Exact jump matrices for P{w) near w ~ —1. Let J7 be a fixed neighborhood of w = —1. By 
straightforward substitutions, we may assume without loss of generality that within U, the jump contour 
for F(w) consists of the real axis together with four arcs (the lens boundaries) meeting at some real point 
w = ui* G J7nM tending to —1 as criticality is approached; the way is determined as a function of x and t 



will be explained later (the resulting formula being (4.52)). We suppose further that the lens boundaries lie 
along certain curves (also to be specified later) tangent at w = to the straight lines arg(u; — w^,) — ±7r/3 
and aig{w — w*) = ±27r/3. Some calculations show that the exact jump conditions for F{w) within U can 
be expressed only in terms of the analytic function T]\[{w) « 1 and another analytic function k{w), which is 
the analytic continuation from w < min{w*, —1} of the function 



k{w) — k(w; x, t) := 2iQ{w; x, t) + L{w) — 2g(w7; x,t), w < —1. 
The jump matrix Vp(^) for which along each of six contour arcs meeting at w = 



(4.1) 

we have P+(C) = 

P_ (C)Vp(^) is illustrated in Figure 4.1 The function Tpf{'w) can easily be removed from the jump conditions 




Figure 4.1. The jump matrix for P(w) near w = ^ —1. The four non-real arcs of the 
jump contour are in general only approximately linear near the intersection point w = w^.. 
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by making the following near-identity transformation (locally, in U) 

P{w) := P{w)dN{wy' 
where (In^w) is the piecewise-analytic function given by 



in the region tangent to the sector | arg(w* — < ^ 
in the region tangent to the sector | arg(w7 — w^)\ < ^ 
otherwise. 



(4.2) 



(4.3) 



Note that djv(w) = 1 + 0{€n) holds uniformly for w £ U, so indeed P{w) ^P{'w) = I + ©(eat). The jump 



conditions satisfied by P{w) near w 
all cases. 



1 are as shown in Figure 4.1 but with Tn{w) replaced by 1 in 



4.2. Expansion of k{w) about w — —1. 

Proposition 4.1. The function k{w) — k{w]x,t) is analytic at w 



-I, having the Taylor expansion 



kiw] X, t) = -t+ -l\x{w + 1) + 



1 A 1 

-Ax 1 

4 8 



[w + lf 



Ax 
Ax ■ 



-t + v 



15 



{w + lf 



(4.4) 
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where v > is independent of w, x, and t, being defined by (1.30). All Taylor coefficients are linear in Ax 
and t, and hence in particular the error term is uniform for bounded Ax and t. 

Proof. For w < min{w*, —1} and w G [/, we have the formula 

1 



k{w;x,t) — 2iQ{w;x,t) 



where m{w) is given by 



m{w) 



t 



[m+{w) +m-{w)], 



%{s) ds 



(4.5) 



(4.6) 



2 2ttJ„ (-s)i/2(s_iy)- 

To express the boundary values m±(w) for w < min{u;*, —1} in a form admitting analytic continuation to a 
full neighborhood of it; = —1 (a point that lies on the discontinuity contour [a, b] of m(u')), let C± denote two 
contours from s = a to s = b, with C-^- in the upper half s-plane and C_ in the lower half s-plane, and such 
that the two contours are mapped onto each other with reversal of orientation under the mapping s i-^ s~^ 
(in particular, this mapping permutes the points o and b). Then, also recalling the definition of Q(w;x,t) 
in terms of E{w) and D{w), 



Kw) = {-w) 



-,„^l/2_(_^)-l/2^ _ 



{-wf'^ + {-w)-''^ 



-w) 



1/2 



2tt 



%{s) ds 



C,UC_ (-.)V2(.-^)- 



(4.7) 



This formula represents the analytic continuation to a full neighborhood of w 



l of the function k 



originally defined by (4.1) only for w < Taiii{w^, —1}. In fact, the domain of analyticity for k has now been 
extended to the domain enclosed by the contour C+ U C_ . It is also obvious that k depends linearly on x 
and t, and this property will clearly be inherited by all of its Taylor coefficients. 

To analyze k{w) near the point w = —1 to which the self-intersection point w^ will converge at criticality, 
we begin with the following elementary Taylor expansions (convergent for jw + 1| < 1): 



i-wf/^ - (-«;)-i/2 
i((-i.)V2 + (-^)-V2) 



liw + l) + \iw + lf + ^iw + lf + ^{w + If + 0{{w + If) 



-I 



1 



1 



{w + lY~-{w + l) 
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[w + lf + 0{{w + lf). 



(4.8) 



We also have the expansion 



s — w 





1 
1 


1 1 

1 1 




iw 


S + i 


[s + ij 




1 


+ 


1 


1 

L 


1 






.(.s + 1 


)3 2{s 


+ 1)2 


+ 


1 


1 


1 






.(.s + 1 


)4 2 is 


+ 1)3 


+ 


1 


1 


1 






[(.s + 1 


)5 2 is 


+ 1)4 


+ 0{{w + 









1 1 



8.S + 1 
1 1 



iw + If 
1 1 



iw + lf 



(4.9) 



2(s + l)4 8(s + l)3 16(s + l)2 128S + 1 



iw + lf 



which is a convergent power series for |w + 1| < min{l, \s + 1|}. Therefore, given contours C± as above, if 
jiy + 1| is sufficiently small we may integrate term- by-term to obtain 

(-w;)i/2 



27r 



9ois) ds 



(w + 


1) + 


h- 






1 












16 






1 












16 



1 



iw + lf 



h h 
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where 



h. ■■= 



2ti 



%is) ds 



c+uC- 



(-S)l/2(s+ l)fc^ 



By making the substitution a — s ^, we may rewrite these integrals in the form 



h 



1 

'2^ 



%{a) da 



c^uc. (-a)V2(a + l)(a-i + l)'=- 
h 



Therefore, we see immediately that 

-h = 0, 

while by averaging (4.111 and (4.12) we obtain for I2, I3, and I4 the following: 

1 



' ' 47r L ("5)1/2(^ + 1)1 + ^^ 



ds 



and 



h = 



1 



1 - s3 



' 4^yc,uc_ (-s)V2(. + l)(. + l)= 



ds. 



iw + iy 

iw + lf + Oiiw + lf), 

(4.10) 

(4.11) 

(4.12) 
(4.13) 
(4.14) 

(4.15) 



Also, averaging as given by (4.11 1 and (4.12 ) and comparing with (4.14) and (4.15 1 we obtain the identity 

(4.16) 



2/4 -h^ 2/4 - A 



Now we calculate I2 — I3 exactly. Recalling the definitions (1.10) of the functions i?(-) and £*(•): as well 
as the definition (1.14) of 6*0(5), we see that 

^{E{s)) 



-E{s)E'is) ds. 



(4.17) 



87ryc+uc_ Disr 

Since Dis) can be eliminated in favor of -E(s) by the identity D{sf = Eisf + j, we would like to introduce 
V — --AiEis) as a new integration variable. The two contours C± are mapped under — 4zi?(-) to two 
oppositely-oriented copies of the same contour. In particular, —AiEiC^) is a teardrop-shaped contour 
beginning and ending at the point v = — 4i£'(a) — — 4i£'(fa) lying on the real axis to the right of w 2 
and encircling the point v — 2 — — 4«i?(— 1) once in the positive sense. Along the contours C± we have the 
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identities D{s) = ±(i?(s)^ + |)^/^ (principal branch). These considerations show that both contributions 
from C+ and from C_ are equal, and so 



h = h = 



27r 



\E'(w/4) 



-{-2v)dv. 



Next, we integrate by parts, using the fact that '^{E{a)) — 0: 



U = - 



ip{v) dv 

^J-MEiC^) (4-«2)l/2^ 



(p{v) 



dv 



*(ii;/4). 



Collapsing the contour to the top and bottom of the interval (2, —4iE{a)) and noting that — 4i_E(o) 
we obtain 



/2 = /3 = 



ipjv) dv 



(4.18) 

(4.19) 

-G(0) 

(4.20) 



Now substituting from the definition (1.8 1 and exchanging the order of integration (see [7| Proposition 1.1]) 
leads to the identity 

h = h = -^0-^-2) = -^Xcrif (4.21) 

Next we will show that 

(4.22) 



l4 = - ^a;crit 



where > is defined by (1.301. Using (4.21) together with (4.14 )-(4.15 1 and the definitions (1.10) of E{-) 
and D(-) gives 



/4+^X„.it=/4-/2=^ 



Introducing v 



ic+uc- D{s) 
-AiE{s) as a new integration variable, as above, yields 



1 

2 "^crit 



1 

27ri 



HE 



*(zi;/4) 



{-2v)dv, 



and then integrating by parts. 



1 1 /" fi"^) 

Note that by the subsitution m — G(s)^, V'(''^) rnay be written in the form 



^(m) 



dm 



y/m - i;2^G(0)2 - m 



< w < -G(0), 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



where the function 5^ is defined in terms of the initial condition G(-) by ( 1.6 ). We need to write this formula 

those -i 

z - G{Qy 



in a way that admits analytic continuation to complex v such as those v e —AiE{C-). To this end, let 

^2 \ -1/2 



T{m,z) :=z(G(0)2-m)-i (^1 
where the principal branch is meant, and then note that 



m - G(0)^ 



^(to) 



dm, 



(4.27) 



(4.28) 



where O is a loop contour surrounding, in the negative (clockwise) sense, the straight-line branch cut of 
T{m,v'^) (viewed as a func tion of m) connecting with G(0)^. In writing the formula (4.28) we therefore 
are using Assumption 1.4 to guarantee analyticity of J# in a neig hborhood of m = G(0)^. With the loop 
contour O fixed, it is clear that (p{v) as given by (4.28) is analytic ii z — v'^ varies in the region enclosed 
by O. In particular, we will assume that G_ has been chosen so that — 4ii?(G_) is completely contained in 
this region. Since —AiE{C-) and L are both fixed contours, we may exchange the order of integration upon 
substituting (|4.28| into (|4.25|: 



I4 + ^■^crit 



247ri 



^(m) 



j(m) dm, 



(4.29) 
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where the inner integral is now: 
j{m) 



T{m, 



„2)3/2 



{2v) dv 



T(m, z) dz 

-4i£;(c_)]2 (4 



7)3/2 



(4.30) 



-itE{C-) (4 ■ 

Note that the contour [— 4ii?(C_)]^ is a teardrop-shaped contour beginning and ending at z = G(0)^ and 
encircHng the part of the (principal) branch cut of (4 — z)^/^ between z — A and z — G(0)^ > 4 once in 
the positive sense, and each m £ O lies outside of this closed contour. To evaluate j{m) we note that the 
branch cut of T(m, z) viewed as a function of z is the ray from z — m to z = 00 in the direction away 
from z = G(0)^. By simple contour deformations of [— 4i_E(C'_)]^ taking into account that the integrand 
T{m, z)(4 — z)~'^/^ is integrable at 2; = 00 and changes sign across all branch cuts, we obtain, for 5{to} 7^ 0, 

/ X r IX Z"™ T(m,z)dz 

jm=2sgn3{™} / ' , 4.31 

Jcioy (4-z)3/2 

where the path of integration is a straight line in the region of analyticity of the integrand. Since sgn(S{m}) 
and sgn(3{2:}) coincide, we may absorb the sign by changing the branch of (4 — z)^/^ in the lower-half 
z-plane: 



G(0)- 



T{m, z) dz 
i{z-Af/^' 



With the substitution 

the formula for j{m) becomes 
j{m) = 



a 



m 



1/2 



(m 



da 



- {m-G{Of)a 
a{m — 4 — (m 



< a < 1, 



-3/2 da 



(4.32) 



(4.33) 



G(0)2)a2)-i/2 



da 



(4.34) 



(4-m)^G(0)2 -4' 

Inserting this formula into (4.29) and taking into account that the only singularity of the integrand enclosed 

4 yields 



by the contour O is the simple pole at m 

h 



1 



^(4) 



(4.35) 



Evaluating ^^(4) directly using the definition (1.6) then proves that indeed /4 = ^ where v > Q is 



given by (1.30) 



sions (4.^ 



Using (4.13), (4.16), (4.21), and (4.22 



in (4.10), and combining the resulting expansion with the expan- 

□ 



and the analytic formula (4.7) for k(w,x,t) completes the proof of the Proposition. 



4.3. Conformal mapping near w ^ —1 and new spacetime coordinates. From (4.4) it is clear that 
exactly at criticality, /c(w) = i'{w + l)^ + 0{{w + l)'^). For small |Ax| and |t|, the cubic degeneration unfolds, 
with one real and two complex conjugate roots, or three real roots, near w = —1. The double critical point 
(double root of fc'(w)) unfolds generically to a pair w = w±{x,t) of simple critical points near w — —1, 
and either these are both real or they form a complex-conjugate pair (in both cases since k(w) =^ k{w*)*, 
{k{w+) - k{w-)f e E). 

The main idea here is that to unfold the cubic degeneracy it is really only necessary to take into account 
the lower-order terms in the Taylor expansion, and therefore it seems desirable to somehow replace k{w,x,t) 
by an appropriate cubic polynomial with coefficients depending on {x,t). This issue has arisen frequently 
in the construction of local parametrices for matrix Riemann-Hilbert problems corresponding to certain 
double- scaling limits. Perhaps the first time such a replacement was made rigorous was in the paper of 
Baik, Deift, and Johansson [2 , in which a certain double-scaling limit of orthogonal polynomials on the unit 
circle is analyzed, and the authors construct a local parametrix by essentially truncating an analogue of the 
Taylor expansion (4.4) after the cubic term. Of course such a truncation is not exact, so there are errors 



incurred in modeling the jump matrices by others having cubic exponents, and the effect of these errors must 
be carefully controlled. Also, any artificial truncation of a Taylor series can only be accurate if the local 
parametrix is constructed in a disk centered at the expansion point whose radius tends to zero at some rate 
tied to the large parameter in the problem. This fact further implies that estimates must be supplied for 
norms of singular integral operators that are independent of the moving contour and it also typically means 
that divergence of the outer parametrix near the expansion point can lead to a sub-optimal mismatch with 
the local parametrix, leading in turn to sub-optimal estimates of errors. A significant advance was made 
by Claeys and Kuijlaars in ■12] (see in particular section 5.6 of that paper); here the authors consider a 
similar double-scaling limit and take the approach of constructing a certain nontrivial conformal mapping 
W — W{w) of a fixed neighborhood of the expansion point to a neighborhood of the origin. The conformal 
mapping is more-or-less explicit, and it does not depend on the parameters (analogues of x and t) driving 
the system to criticality. However, the point is that in a full fixed neighborhood of the expansion point, 
the analogue of the function k{w]x,t) is represented exactly (no truncation required) as a cubic polynomial 
in W . The difficulty that remains with this approach is that the coefficients of the cubic in W actually 
depend on w (so in fact it is not really a cubic after all). The remarkable approach of [T2] is to simply 
substitute the w-dependent coefficients into a known solution formula for exactly cubic exponents (involving 
solutions of the linear differential equation whose isomonodromy deformations are governed by solutions 
of the Painleve-II equation), resulting in a local parametrix in terms of Painleve transcendents depending 
(through the coefficients of the cubic) on w. Confirming that such an approach actually provides a usable 
local parametrix requires exploiting a priori knowledge of the behavior of solutions of the nonlinear Painleve- 
II equation and its auxiliary linear differential equation. If this information is available in a convenient form, 
the approach of Claeys and Kuijlaars delivers a vast improvement over earlier methods because it really works 
in a neighborhood of fixed size centered at the expansion point where critical points coalesce at criticality. 
The technique advanced in |12j has more recently been applied to problems of nonlinear wave theory as well 
(see, for example, [TO]). 

We choose instead to replace k{w,x,t) by a cubic in an exact way, an approach that provides all of 
the accuracy of the Claeys-Kuijlaars method but seems simpler and requires no a priori knowledge of 
the behavior of solutions of the local parametrix Riemann-Hilbert problem. The approach we are going 
to explain now has also been used recently to study a different kind of double-scaling limit for a matrix 
Riemann-Hilbert problem in [5]. As part of a careful study of the asymptotic behavior of exponential 
integrals with exponent functions having coalescing critical points (to generalize the steepest descent or 
saddle point method), Chester, Friedman, and Ursell [S] showed how to construct a substitution that rendered 
the exponent function in the integrand in the exact form of a cubic polynomial. Their method applies in 



the current context to establish that, because the coefficient v in the expansion (4.4) of k{w; x,t) is strictly 
positive, there is a suitable choice of new spacetime coordinates r — r{x,t) and s — s{x,t) for which the 
relation 

k{w;x,t) ^W^ + rW - s (4.36) 

defines an invertible conformal mapping W = W{w) = W{w; x, t) between w G U and W G W{U) that 
preserves the real axis, for (x, t) near enough to criticality. Moreover, the new coordinates r and s depend 
continuously on {x,t) near criticality. Unlike in the Claeys-Kuijlaars method |12| . neither r nor s depends 
on w with the cost that the conformal mapping will now depend on {x,t). 



The new coordinates r — r{x,t) and s — s{x,t) are to be determined so that under (4.36) the critical 



points of the cubic on the right-hand side, namely W = ±(— r/3)^/^, correspond to the two critical points 



of k{w) near w — —1 when Ax and t are sufficiently small. Evaluating (4.36) for w ~ 'w±{x,t) and 
W = ±{—r{x,t)/3y^^, one obtains the formulae 



r{x,t) = -3 [ — [k{w+{x,t);x,t) - k{w-{x,t); x,t)f 
s{x, ^) — ~2 2;, t) + k{w-{x, t);x, t)] , 



(4.37) 



where in the formula for r the real cube root is meant, and hence both r and s are real. Moreover, it is 
possible to show that r(x, t) and s{x, t) are analytic functions of x and t near criticality, and have two- variable 
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Taylor expansions of the form: 

r{x,t) = -\j-Ax + 0{Ax^,tAx,f) 
^ ' 2i^i/3 ^ ' ' ' (4.38) 

s{x, t) = t + 0{Ax^,tAx, t^). 

To see this, first we find a unique root w = —1 + a of k" {w; x, t) such that a = at criticality and such that 
a — a(x,t) is an analytic function of (x,i) near criticality. Indeed, from the Taylor expansion (4.4) and the 
fact that v ^ we see that the analytic implicit function theorem applies and we obtain 

a(x, t) = — I- Ax + J-t + 0(Ax'^,tAx, f). (4.39) 



Now write w = —1 + a{x,t) + z, and re-expand k(w;x,t) about 2 = 0. From (4.4 1 and the definition of 
a{x,t) we obtain the convergent power series expansion (note that there is no quadratic term) 

oo 

k{w]x,t) = fco(a;,t) + ki{x,t)z + ^ fc„(a;,t)z", (4.40) 

n— 3 

where the coefficients are all analytic functions of (x, t) near criticality, and in particular, 

ko{x,t) ^ -t + 0{Ax^,tAx,f), (4.41) 

ki{x, t) = ^Ax + 0{Ax^,tAx, t^), (4.42) 

k3{x,t) = iy + 0{Ax,t). (4.43) 

The radius of convergence of this series is bounded away from zero near criticality. The equation satisfied 
by the critical points of k is 

CO 

k'{w;x,t) = ki{x,t) + 3k3{x,t)z'^ + ^nkri{x,t)z'''^ =0. (4.44) 

ri=4 

Let p be any number satisfying the equation 

P" = -¥r^V (4.45) 
3k3[x,t) 

Note that is an analytic function of {x,t) near criticality because fc3(xcrit,0) = > 0. If p = then 
ki{x,t) = and z = is a double root of ( |4.44 ). So we suppose that p 7^ 0, and we rescale z by writing 
z = pq for some new unknown q. Dividing through by 3fc3(x, t) ^ and canceling a factor of p^ then converts 
(4.44) into the equation 

+ E ^^^P"-'?"-^ = 1- (4-46) 

The coefficients of g"^^ in the sum are analytic functions of the three variables x, t, and p. When p = 
there are two solutions, q — ±1, and by the implicit function theorem there are two corresponding solutions 
for p w and {x,t) near criticality. The two solutions are related by the symmetry {p,q) 1— > {^p,—q). We 
may write them in the form 

9 = g± (x, t, p) = ± + g„ (x, t) (±p)" j . (4.47) 

The coefficients qn{x,t) are all analytic functions of {x,t) near criticality, and the radius of convergence of 
the series with respect to p is bounded away from zero near criticality. The corresponding critical points are 
written in terms of z as 

00 

z = z±{x, t,p)=±p + Y. *)(±/')"^'- (4-48) 

n=l 
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To calculate r{x,t) and s{x,t) we need to evaluate k{w,x,t) at the critical points. This is accomplished by 
substitution of the series (4.48) into (4.40): 

k{-l + a{x, t) + z± {x, t, p); X, t) = fco(a;, t) + fci(x, t) {±p + (7„(x, t){±pY+^^ 

oo / oo \ ™ 

+ k^{x, t)i±p+J2 Inix, t)(±pr+^ 

m—3 \ 71—1 J 

oo 

= fco(a:,i)±p^^5n(a;,i)(±Pr, 



(4.49) 



Tl = 



where we have used (4.45 ), and 50(2^1 i) •= ^2fc3(a;, i). Here again, the coefficients gn{x, t) are certain analytic 
functions of (x^t) near criticality, and the radius of convergence of the power series in ±p is bounded below 
near criticality. Then, by definition, we have 

2\ 1/3 



ki{x,t) 



2g^{x,t)p^ [l + Y. 



\ goix,t) 



k3{x,ty/^ 



1 - 



^ 00 



92j{x,t) 
{ kz{x,t) 



ki{x,t) 
Sks^x, t) 



2/3 



(4.50) 



where we have used the definition of go{x,t) in terms of k^^^x^i) and ( 4.45[ ). Therefore, r{x,t) is clearly an 
analytic function of (x, t) near criticality. Similarly, by definition we have 



s{x, t) = — 



1 



2fco(.x,t) + 2£52,+i(a:,t)(±p)2^'+4 



= -ko{x,t) - ^52j+i(a;,t) I - 



ki(x,t) 
3k3{x,t) 



3+2 



(4.51) 



which again is obviously an analytic function of (x, t) near criticality. The leading terms of r(x, t) and s{x, t) 
near criticality are easy to calculate from these formulae, with the result being (4.381. 

As W{-; X, t) is a conformal map, its inverse is an analytic function mapping a neighborhood of = to 
a neighborhood of m; = —1. We now define 

w^^w^{x,t)::=:W-^{Q;x,t), (4.52) 

which is a real analytic function of x and t near criticality. Note also that W' {w■^,{x,t)]x,t) > (prime 
denotes differentiation with respect to w) near criticality. It is not difficult to obtain the Taylor expansion 
of the conformal mapping W{w;x,t) about w — exactly at criticality. Indeed, from (4.4) and ( |4.36 1 with 
Ax = t — r — s — 0, we obtain the equation 



l + -(u; + l)+0((w + l)2) 



and analytically blowing up the cubic degeneracy we obtain 



W{w) = t/i/3(w + 1) 
In particular, this implies that 



l + ^(w + l)+0((«; + l)2) 



at criticality. 



, at criticality. 



(4.53) 



(4.54) 



W'{w^) = W"{w^) = v^^^ > 0, at criticality. (4.55) 

Also, it is clear that 1; Xdit: 0) = 0, and therefore w*(a;critjO) = —I. Since w^,{x,t) is an analytic 
function of x and <, it then follows that 



w^xA) ^ -l + 0{Ax,t), 
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(4.56) 



and furthermore, 

W'{w,{x,t);x,t) = i^^/^ + 0{Ax,t) and W" {w,{x,t); x,t) = i^^/^ + 0{Ax,t) (4.57) 
near criticality. 

Now let ( and y be scaled versions of W and r respectively: 

W r 

We will regard y as being bounded. It then follows that without approximation the exponent appearing in 
the jump matrix 'Vp{w) takes the form of a simple cubic with a formally large constant term: 

'^^^C^+.C--. (4.59) 

4.4. Formulation of an inner model problem. To begin with, we wish to find a simpler representation 
for P™*(-u;) valid when w is close to w = sa — 1. Using (4.58), we write P™*(w) in the form 

P^H = Vm{wr{-W{w))("'~'^'^^- = r;„(i.)-^'e(^"-i)'^^/^(-C)(2"-i)-^/2, (4.60) 

where rim{w) is independent of e^r and is analytic and nonvanishing in a neighborhood U oi w ^ w■^,^. 

Tjr/ \ \ 1/2— m 

W{w) ^ ' 



iV^ + V-w.)'-'"' • (4.61) 

W — W^. 



(Analyticity follows since from (4.52) we have W{w^) — and W is analytic at with W' [w^) > 0.) 

Given a value of to G Z (to be determined below), we wish to construct an inner model, valid for w €z U, 
of the matrix P{w) (related to P{w) ioi w G U via (4.2)). We temporarily denote this model as W(C(ti;)), 
and we want it to have the following properties: 

• Supposing that in U the lens boundaries are identified with arcs of the curves a,Tg{W (w)) = ±7r/3 
and a.Tg{W {w)) — ±27r/3, which makes them segments (of length proportional to ) of straight 
rays in the ^-plane, the inner model W(^(ii;)) should be analytic exactly where F{w) is within U 
and should satisfy exactly the same jump conditions that P{w) does within U. 

• The inner model W(C) should match onto the latter factors in (4.601 along the disc boundary dU in 

1/3 

the sense that W(C) may be analytically continued from each sector of the domain C G e^r W{U) 
to the corresponding infinite sector in the C-plane, and that 

lim W(C)ei^^""^"^/'(-C)^'"'™'"^/' = I, (4.62) 

with the limit being uniform with respect to direction in each of the six sectors of analyticity. Since 
dU is fixed and bounded away from « —1, upon scaling its image under W{-;x,t) by to 
work in terms of the variable ^, we see that w G dU corresponds to C — > oo at a uniform rate of 

Since W(C(w)) depends on w through it is convenient to use (4.59) to write the jump matrices for 

F{w) (and hence also for W(C(i«))) in terms of which shows that the jump matrices involve the product of 
exponentials e^*/'^"e*('> +2''^). The constant (^-independent) factors e^*/^" present in the jump matrices can 
be removed, and simultaneously the presence of ejv in the normalization condition (4.62) can be eliminated, 
by defining the equivalent unknown 

Z,„(C;j/) := 4-2m)a3/6^_|..3/,„^^^^gi,^3/e„_ ^4_g3^ 

The conditions previously discussed as being desirable for W(C) are then easily seen to be equivalent to 
those of the following problem for Zm{C, y)- 

Riemann-Hilbert Problem 4.1 (Inner model problem near criticality). Let a real number y and an integer 

TO be fixed. Seek a matrix Zm(Cj 2/) with the following properties: 

Analyticity: ZmiCw) analytic in C except along the rays arg(C) = kn/S, k — 0, ...,5, from 
each sector of analyticity it may be continued to a slightly larger sector, and in each sector is Holder 
continuous up to the boundary in a neighborhood of Q — 0. 

28 



Jump condition: The jump conditions satisfied by the matrix function Z„i(C; y) are of the form 
Zm+{C',y) = Z„i_(C; ?/)Vz(C; y)j where the jump matrix "V y) is as shown in Figure ^.2 and all 




Figure 4.2. The jump matrix Yz{C;y)- 



rays are oriented toward infinity. 

Normalization: The matrix 7i„i{C.; y) satisifies the condition 

lim Z„(C; y)(-C)''"'™'"^/' = I (4.64) 
with the limit being uniform with respect to direction in each of the six sectors of analyticity. 
Suppose that Riemann-Hilbert Problem 



4.1 



has a unique solution for y e {y-,y^) C K and for some 
integer to. We will now describe how to use it to create a model for P(w) valid in the neighborhood U 
assuming that y e {y^,y^); we set 



(4.65) 



The effect of post-multiplication by dNlw)^"^^ is simply to restore the factors involving Tn{w) to the jump 
conditions. We emphasize that within the fixed neighborhood U of w = — 1, the inner model matrix P"(?i') 
satisfies exactly the same jump conditions along the six contour arcs meeting a,t w — as does F(w). It 
appears reasonable to propose a global model for P{w) in the following form: 



w eU 



(4.66) 



The integer m is evidently at our disposal. We will later describe how it should be chosen. 
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5. Solution of the Inner Model Problem 



5.1. Symmetry analysis of the inner model problem. While it may not yet be clear whether there 
exists a solution of Riemann-Hilbert Problem |4.1| for any y G M and m € Z at all, it is a standard argument 
that there exists at most one solution for each y G M and m G Z, and that every solution must satisfy 
det(Z,„(C; y)) = 1. In this short section we suppose that y G M and m G Z are values for which Riemann- 
Hilbert Problem 4.1 has a (unique and unimodular) solution, and we examine some of the consequences 
of the existence. We will later show that Zm(C; 2/)(~C)^^~^™'''^^^^ has a complete asymptotic expansion in 
descending powers of C as C — oo: 



(5.1) 



and where the matrix coefhcients are the same regardless of which of the six sectors of analyticity of Z,„(C; y) 
is used to compute the expansion. 

Given a solution Zm(C; y) of Riemann-Hilbert Problem 



4.1 



consider the related matrix 



F(C) :=Z™(C;y)* 



(5.2) 



Since the jump contour for Z„i(C;y) is invariant under complex conjugation, F((^) will be analytic in the 
same domain that Z,„(<^; y) is. It is a direct calculation that F((') satisfies exactly the same jump conditions 
as does 'Lm{C,] y)- Therefore, the matrix F(C)Zm(<^; y)~^ extends continuously to the jump contour from both 
sides of each arc. It follows from the classical sense in which the boundary values are attained (even at the 
self-intersection point C — 0), and from the fact that Zm(C;y) is unimodular, that in fact F(C)Zm(C;y)^^ is 
an entire function of C- Moreover, 



lim F(C)Z„(C;y)- 



lim F(C)(-C)^^"^"^'''/^ 



z,„(C;y)(-C)^'-'"^"-'/' 



= hm Z,„(C;y)*(-C)('-'")"-^/'= Hm Z„,(C*; y)(-C)('-'")"^/' 



(5.3) 



so by Liouville's Theorem, Zm(C*;y)* = F(C) = Zm(C;y). Applying this symmetry to the expansion (5.1) 
shows that the elements of the matrices Am(y) and 'Qmiy) (and in fact all of the matrix expansion coefficients) 
are real numbers. 

We can obtain more detailed information by exploiting a further symmetry in the special case of to = 0. 
Indeed, let us compare Zo(C;y) with the matrix 



F(C) := 



Zo(-C;y)*cri, 
-Zo(-C;y)icri, 



3{C} > 

3{C} < 0. 



(5.4) 



Clearly, det(F(C)) = 1, and it is a direct matter to check that F(C) is analytic precisely where Zo(C;y) is, 
and F(C) satisfies exactly the same jump conditions on the six rays of Riemann-Hilbert Problem |4.1| as does 
Zo(C;?/)- Therefore, F(^)Zo(C; y)~^ is an entire unimodular matrix function of C,. To identify this entire 
function as a polynomial, it is enough to extract the non-decaying terms in its asymptotic expansion as 
C — oo; Assuming without loss of generality that 3{C} > so that (— C)"^'^^^ = ^^o'aC'^^^^j we have 



F(C)Zo(C;y)-^=F(C)(-za3)r^/' 

= Zo(-C;y)^7i^T3r'/' 

= Zo(-C;y)r^/' 

= (l-Ao(y)r^ 
= CT-C- Ao(y)cr- 



-c 



-A,{y)C'+0{C^)) 

(I-Ao(y)r^+O(r')) 



(5.5) 



-Q-^' 
C 

-a_Ao(y) + 0(r^), C 



(I-Ao(y)r'+O(r')) 



where 



1 




and 
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1 



(5.6) 



and so F(C)Zo(C;2/) ^ may be identified with the polynomial terms in this expansion; therefore we arrive at 
the identity 

Since both sides represent unimodular matrices, we learn that 

AoMvf - 1- (5.8) 



We will later obtain a refinement of (5.8), namely that ^0^12(2/) = 1 holds for all real y. The sign information 
will follow from the explicit solution of Riemann-Hilbcrt Problem 4.1 for m — in terms of special functions. 



which will be possible precisely as a consequence of the identity (5.8 1 



5.2. Differential equations derived from the inner model problem. Lax equations and Painleve- 

II. The only parameters in the Riemann-Hilbert problem characterizing Zm(C;?/) are m e Z and y € 
M. Assuming existence for some fixed m € Z and for y in some open set, we will now investigate some 
consequences of the dependence of the solution on y. We suppose that the expansion (5.1) holds in the 
stronger sense that the following are also true: 



d_ 
dy 



Z™(C;2/)(-C) 



_A\(l-2m)cr3/2 



= 0{C^) and 



d_ 



Z™(C;2/)(-C) 



(l-2m)o-3/2 



= 0(0, c 



(5.9) 



The matrix Lm(y, defined by 

L™(y,C):=Z™(C;2/)e-('^'+^^)'^^/2 (5.10) 
has jump discontinuities along the six rays mediated by jump matrices that are independent of both y and 
C,. Given the compatibility of the jump matrices at C = 0, this implies that the matrices 



U™(2/,C) := ^^Lj and Y^{y, () ^^Kn 



(5.11) 



are both entire functions of C- Using the above expansions of Zm(C;j/) and its derivatives, we easily obtain 
the following expansions involving lim- 



9L„, 
dy 



9L„, 



L-1 



de 

dy 
1 



-0(73 



1 



= -^CT,C--A,r,iy)a3 + 0{C') e 



2 

de 



,-6(73 



[I + A„,{y)C^ + B„,{y)C^ + 0{C'')] + 0(C'') e 



,-0(73 



(5.12) 



= -o'^sC' - ^A„,(y)a3C - ^B„,{y)a3 - -ya^ + 0{C') e 



,-0(73 



„0(73 



where 



e:=^[C'+2/C+(l-2m)log(-C)] 



(5.13) 



We can then easily see that both Um and Vm grow algebraically in ^ as C — ^ c)o and hence are necessarily 
polynomials: 



U™(y,C) 
V™(y,C) 



1 



3 3 3 3 

'2<^3C^ + -[(73, A„(y)]C - -[<73,A,niy)]Amiy) + -[(73,Brniy)] - ^yCTs- 



(5.14) 



We may rewrite these in the form 



U,„(y,C) = -^^3C + 



V™(y,C) = -2'T3C' + 3 



U^{y) 

-Vrniy) 

U,n{y) 
-Vm{y) 



c + 



-6Urn{y)Vrn{y) - y 2W„(y) 

-2Z„,(y) 6Umiy)Vmiy) + y 



(5.15) 



where 



V,„(2/) 
W,„(2/) 

Zrn{y) 



A. 



m,21 



iy) 



= iBm^iiy) ~ iArn,2i{y) Am, 11 (y)- 



The matrix L,„ is therefore a simultaneous (and fundamental, since det(Lfrj 
equations 

U™(y, C)L™ and = V™(?/, C)L„, 



9y 

and therefore the compatibility condition 



dV 

9y 



5C 



[U„,,V,„] =0 



(5.16) 

1) solution of the Lax 
(5.17) 

(5.18) 



holds identically in C,. Separating the powers of C, leads to the following system of differential equations (the 
Painleve-II system) governing the quantities defined by ( |5.16 1: 

^Liy) = -l^rniy) 



vLiy) = ^Zmiy) 

W'miy) = 6Ur,,{yfVm{y) + yU,n{y) 
- -(sUm{y)Vm{yf - yV,„(y). 



(5.19) 



Although satisfied by quantities evidently depending on m e Z, the Painleve-Il system does not involve m 
in any explicit way. Eliminating Wm and Z^, yields the coupled system of second-order Painleve-II-type 
equations 

U';^{y) + 2Um{yfVm{y) + \yUm{y) = and V:^{y) + 2Um[y)'Vm{yf + \yVm[y) = 0. (5.20) 

In fact, functions related in an elementary way to Um and Vm turn out to satisfy uncoupled second-order 
equations of Painleve-II type, but we must use the inhomogeneous form of the Painleve-II equation, and the 
coefficient of inhomogeneity will depend on m € Z. To see that this is so, we follow [18, pages 154-155] by 
multiplying the first equation of ( |5.20 1 by Vm , the second by Um , and subtracting: 

= Vm{ymy)-Um{y)v:'M = ^ [Vm{y)U'm{y) -Um{y)Vm{y)] ■ (5.21) 

ay 



Hence, the quantity 



1 



A,„ := Vm{y)U'miy)~i^m{y)v:jy) = -- [Vm{y)Wm{y) +Um{y)Zm{y)] 



(5.22) 



is a constant, independent of y. 

To compute the value of Am for the particular solution of (5.19) corresponding to the specific Riemann- 
Hilbert problem for Zm(C;y), we first consider the direct problem for the Lax equation lim,c — V„jLj„, 



formally expanding Lm(y, C) ^r large C by assuming the form 



Lm(2/,C) 




Ay)C 



c 



(5.23) 



where each of the coefficients Y„ i.(?/) is an off-diagonal matrix, and where Dm(2/, C) is a diagonal matrix 
with the asymptotic expansion 



D„(y,C)~D„,_3(y)C' 



_2(y)C'+Dm,_l(2/)C + D™,o(y)log(-C) + ^D„,fe(^/)r^ C^oo, (5.24) 



k=l 
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where all of the coefficients 'Dm,k{y) are diagonal matrices. The value of Am will be determined from the 
coefficient Dm,o- By differentiation of the formal series, we obtain 



dC 



where the matrix coefficients are systematically determined as follows: 



M„ _2 3D,„ _3(j/) + -(Tg 



(5.25) 



3Y„,i(y)D™,_3(y) + 2D™,_2(j/) - 3 



U,n{y) 
-Vmiy) 

M,„,o := 3Y„^2(y)D,„,_3(y) + 2Y,„,i(2/)D„,_2(y) + D,„,_i(y) 



:CT3Y™,i(y) 



-IK^,„(y)V„(y)-j/ 2W„(y) 

-2Zm{y) &Um{y)Vm{y) + y 



- 3 



U,^{y) 
-V,n{y) 



■■■ m.l {y) + 2^3Y™,2(y) 



M„,i := 3Y„,3(y)D™,_3(y) + 2Y„,2(2/)D„,_2(y) + Y™a(y)D™,_i(2/) + D™,o(y) 



-(K^„(y)V,„(y)-y 2W„(y) 

-2Z„(y) &U^{y)V^{y) + y 



^rnAy) 



Z^,„(y) 
-Vmiy) 

-Dm,i(y) ~ ^ 



Ym,2(y) + 2^3Y™,3(y) 



2W,„(y) 
-2Zraiy) 



Z^™(y)' 
-V™(2/) 



Y™.2(y) 

Ym,3(y) + off-diagonal terms, 



(5.26) 



and so on. By setting each 'M.m,k to zero in sequence, and by separating the diagonal and off-diagonal parts 
of each of these equations, we obtain the following: the matrix equation M„i__2 = implies 



D 



m, — 3 



(y) 



1 



70-3, 



the matrix equation _i = implies 

Ym,i(y) = 

the matrix equation Mm.o = implies 



U„,{y) 
V„(y) 



Ym,2(2/) 



1 



Wm{y) 
Zm{y) 



and D™ _2(j/) = 0, 



and Tira.-i{y) ^ ^2^^^^' 



(5.27) 



(5.28) 



(5.29) 



the matrix equation M„j.i — implies 
Y,„,3(2/) = - 



^ra{y?Vra{y) + yUra{y) 

3U^{y)Vrniy? + yVraiy) 



and T>m,o{y) = -3A,„cr3, (5.30) 



where we have recalled the definition (5.22) of A^, and finally, the diagonal terms of the matrix equation 
Mm 2=0 imply that 



Dm,i(2/) = -2Hm{y)(J3 



where the Painleve-ll Hamiltonian is 



Hm{y) ^W„,{y)z„^{y) - ^U^{yfv,n{yf - \yUUyWm{y)- 
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(5.31) 
(5.32) 



It follows that the formal asymptotic expansion of L,„(2;, C) for large C takes the form 



Lm(y,C)e 



(C^+yC)<T3/2 



(-0 



3A„CT3 



[Y„,i(2;) + D™a(y)]r' 



-oic 



Comparing with (5.1 1 and (5.10), we see that 



C oo. 

1 m 

^ 6 ^ y' 



and also that 



2 (5.33) 

(5.34) 
(5.35) 
(5.36) 

BmMy) = ^Wm(y) + 2H,n{y)U^iy) and B^^2iiy) = ^Zmiy) - 2i/™(y)V„(y). (5.37) 

Now that the value of Am has been determined according to ( |5.34 1, a direct calculation shows that the 
logarithmic derivatives 

and 



A„(y) = Y„,i(?;) + Dm^i(j/) 

B,„(?/) = Y„^,2{y) + D™,2(y) + ^T>„,A{yf 



-■■ m.l (2/)Dm,i(y)> 



so that in addition to the relations (|5.16|) we have 
and 



1, 



An,iiiy) = -2H„i{y) and ^^,22(2/) = '2Hmiy) 

1 



Um.{y) 



satisfy the uncoupled Painleve-II equations 



v:,,{y) = 2Vm{yf^T;yVm{y) 



1 



V„^{y) 



Ql(y)-2Qm(2;)3 + -yQ„(y) 



(5.38) 



(5.39) 



= 27'„,(2/)3 + -2/7'„(y)--m 



3" 3 

Note that for general m € Z, both equations are of inhomogeneous type 



2Q™(2;)' + -2;Q„,(y) + -(m-l). 



5.3. Solution of Riemann-Hilbert Problem 4.1 for m — Q. Recalling the relation (5.8) and the defini- 
tions (5.16) of the potentials Um, V„i, 'Wm, and Zm, we learn that when m = 0, Uo{y) = cr, where (assuming 
for the moment continuity of Uo{y) with respect to y) cr = ±1 is a fixed sign to be determined. According 
to the system of differential equations (5.19) necessarily satisfied by the potentials, we can determine the 
values of the remaining potentials in the special case of m = 0: 



Uo{y)^<j, Voiy) 



6 



y, Wo(y) = 0, Zoiy) 



a 
2' 



(5.40) 



Therefore, from (5.15) and (5.17) we see that in the special case of m = 0, the linear differential equations 
simultaneously satisfied by the matrix Lo(a, C) — Zo(C; 2/)e~'-^ +yQ<^zl'^ take the form 



dy 



1 



IC 



Lq and 



3 a2 
2^ 



3aC 



(5.41) 



Let 1 = (^1,^2)^ denote either of the columns of Lq. From the first of the above two Lax equations, we 
observe that 



and then that 



1 

6173 



y + 7^C 



(5.42) 
(5.43) 



where F is any solution of Airy's equation F"{^) — £,F{^). Therefore 1 necessarily has the form 

.c7(6-i/3F(0 + K^(0). 
34 



1 



(5.44) 



and by substitution into the second equation of the Lax pair we learn that the Airy function F may only 
depend on C through the variable ^, in other words, we must select for F a linear combination of Ai(^) and 
Bi(^) with coefficients independent of (. 

Now we try to select the correct solutions of Airy's equation to construct Lo(?/, C) and hence Zo(C; y) in 
each of the six sectors of analyticity in such a way that we obtain the desired asymptotics for large C and 
satisfy the required jump conditions. This procedure will also determine the correct value of ti = ±1. 

Let us consider three specific choices of solution to Airy's equation leading to the three solution vectors 



1" 



Ai(C) 



1^ 



Ai(e±2^'/30 
(7(6-i/3e±2"/3Ai'(e±2"/3^) + iCAi(e±2-/3^)) 



[a(6-i/3Ai'(0 + KAi(e)) 
Each of these is an entire function of (. From the standard asymptotic formulae 



(5.45) 



Ai(z) 



-1/4 -2z^''^/3 



+ and Ai'(z) 



,1/4^-22:3/^3 



2^ 



(1 + 0(1^1-3/2)) (546) 



valid as z — i' oo with | arg(z)| < tt — S for any J > 0, we see that (all of the following asymptotic statements 
assume y is held fixed) 



1" 



-l/2„-{C''+yQ/2 



Similarly, 
1± = 



-ere 



±27ri/3 



(48)1/6^ 

(_()i/2e(C'+aC)/2 
^'^ (48)1/6^ 

(g±j7r/3^)l/2g(C='+2/C)/2 

(48)1/6^^ 



o(C-i) 



c 



|arg(C)|<--5, 



C->oo, |arg(-C)|< 



1^ 



(gT2^i/3^)-l/2g-(C+aC)/2 

(48)1/6^ 



o(C-i) ■ 



I arg(e 



±i7r/3 



01 < 



C^oo, |arg(e=F2W3^)|< 



S. 



(5.47) 
(5.48) 

(5.49) 
(5.50) 



We may now define a candidate for Zo(C;2/) that satisfies the required analyticity and normalization prop- 
erties among the conditions of Riemann-Hilbcrt Problem |4.1[ 



Zo(C;y) := 



l/6V7F(^1^^7e-2W3l-)e(C'+yC)^3/2^ 
i/6V^(e"/6l+,(Te-2W3i-)e(C^+yC)-3/2^ 
i/6^(eW6i+^ _^lO)g(C3+yC)-3/2^ 
i/67^(e-"/6l-,-(Tl")e(^'+?'«)"3/2^ 
i/'^0r(e~*"'/6l-,cre2'^^/3l+)e(^'+2"^)'^3/2^ 

l/6^(_il0^ ^g2W3l+)g(C'+yC)^3/2^ 



arg(C)e (0,f) 
arg(C)€ (f,f ) 
arg(C)€ (f ,^) 
arg(C) € (-TT, 3 
arg(C)e(-f ,-f) 



27r) 



(5.51) 



arg(C)€ (-1,0). 

It only remains to verify the jump conditions of Riemann-Hilbert Problem |4.1[ with the help of the identity 

Ai(z) + e-2'^*/3Ai(ze-2'^'/3) + e^'^'^^Aiize^'^'/^) =0, z £ C, (5.52) 

which implies the vector identity 1° + ^ g27ri/3j+ _ q Uging this fact, it is a direct matter to 

confirm that the required jump conditions are indeed satisfied, provided we make the choice of sign a = +1. 
We have therefore proved the following. 



Proposition 5.1. Let m — 0. For every y G M, Riemann-Hilbert Problem \4.1\ has a unique solution given 
explicitly by the formulae (5.45) and (5.511 with a — +1. Moreover, Zo(C;y) is an entire function of y, 



2o(C; ?/)(~C)'^^''^ a complete asymptotic expansion as C —> oo in descending integer powers of C, (as in 



(5.1 ) ) differentiate term-by-term with respect to both y and C,, and the corresponding potentials lAa{y), Vo(y), 
yVo(2/) md Zgly) are given by (5.401 with a = +1. 
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5.4. Schlesinger-Backlund transformations. Solution of Riemann-Hilbert Problem 4.1 for gen- 
eral m e Z. Now we describe an inductive procedure for obtaining Z„j(C; y) uniquely from tlic conditions of 
Riemann-Hilbert Problem |4.1| for m ^ 0. Suppose ?/ € M is a value for which Riemann-Hilbert Problem |4.1| 
has a (unique) solution Z„j(C; y) for some to e Z. Following [THl Chapter 6], we claim that 

Z™±i(C;2/) = (S±,i(2/)C + S±,o(y))Z,„(C;2;) (5.53) 

if the matrices Q(y) and ^{y) can be properly chosen. 

It is immediately clear that since the the prefactor is entire in the right-hand side of (5.53) has the 
same domain of analyticity, achieves its boundary values in the same classical sense, and satisfies exactly 
the same jump conditions, as does Zm(C;y). Of course the analyticity and jump conditions in Riemann- 
Hilbert Problem 



4.1 



are independent of m, so it only remains to impose the normalization condition (4.64) 



on Zm-|-i(C; y) given by (5.53), and this condition takes the form 

hm (S±.i(y)C + S±,o(y))Z™(C;y)(-C) 



;i-2m=F2)cr3/2 



(5.54) 



Assuming the expansion (5.1) (which at the moment is justified only for to = 0), we are requiring that 



(S±,i(y)C + S±,„(z/))(I + A^iy)C' + B^iy)C' + C^{y)C^ + O{C^)){-0 



Tf3 



i + o(C 



Equating the polynomial part of the left-hand side to the identity yields the conditions 



C oo. 
(5.55) 



(S+.i(y)A„,(2/) + S+.o(y)) 



(S+i(y)B„(y) + S+,o(y)A,„(y)) 




1 



S+.i(2/) 





1 

0' 

1 

1 0' 









(5.56) 



which can be solved uniquely for the matrix elements of ^{y) and o(y)- 
S™,i(2/) 



S+,o(2/) = 



-1 


Brn,l2{y) I Am.l2{y) " ^m,22(y) ^m,12(?/) 
-lMm.l2(y) 



W„(y)/(3Z^„(y)) Um{y) 
~l/U^{y) 



(5.57) 



provided that Umiy) 7^ 0, where we have used (5.16), (5.36), and (5.37). Similarly, from (5.55) we obtain the 
conditions 



(S-.i(y)A™(y) 



1 




+ S- i(y) 



1 



1 0' 









(5.58) 




1 



(S-^i(y)B„(y) + S-^o(y)A„(y)) 
which can be solved uniquely for the matrix elements of ^{y) and Q(y) 
S™,i(2/) = 





-1 



S™.o(2/) 



-lM™,2i(2/) 

An,2l{y) Bm,2l{y) / AmMy) ~ 



-1/V™(y) 
V™(y) Z™(y)/(3V™(y)) 



(5.59) 



provided that Vm{y) 7^ 0. When it makes sense, the formula (5.53) subject to (5.57) or (5.59) constitutes a 
so-called discrete isomonodromic Schlesinger transformation. 

The recurrence relation (5.531 also implies a corresponding recurrence for the potentials lArm Vm, Wm, 
and Zjn- Indeed, expanding the right-hand side of (5.551 as I-|-Am±i(y)C^^ +0(C~^) the terms proportional 
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to C "'^ yield the conditions 

(S+_i(2/)C™(2/) + S+o(y)B„(y)) 

and 

(S,-,i(y)C™(t/) + S„^o(y)B™(2;)) 





1 

1 0" 




-(S+i(j;)A™(2/) + S+o(y)) 
(S,-,i(y)A™(y) + S„^o(y)) 



1 



0" 

1 



-A„+i(y) (5.60) 
-Ara-iiy). (5.61) 



Recalling the definition ( 5.16| ) of the potentials {Um, Vm) in terms of the matrix entries of A„ 



the (2, l)-entry of the matrix equation (5.60) and easily obtain the relation V„j_|-i(j/) 
the differential equations (5.19) satisfied by the potentials for all n, we obtain 



we examine 
Then, from 



Um+iiy) 



1 



y^miy) 



1 



which together with the implied relations for Wm+i and Z^+i 



+ -U'^{y) and V™+i(y) 



1 



(5.62) 



Um{y) 

constitutes a Bdcklund transformation for 



the Painleve-II system (5.19). Similarly, the (1, 2)-element of (5.61) together with (5.19) yields the Backlund 
transformation 



L(y) 



1 



Vm{y) 



and V„ 



v{y) = \v:M 



1 



Vm{y) 6 



yV,„(y), 



(5.63) 



which is easily seen to be inverse to (5.62). The iterative scheme for generating the solution to Riemann- 
Hilbert Problem |4.1| for general m G Z is then simply to start with the solution valid for all ?/ S M for 
m = as described by Proposition |5.1| and attempt to apply the discrete isomonodromic Schlesinger 
transformation ( |5.53 ) in order to repeatedly increase or decrease the value of m in integer steps. At 



the level of the implied Backlund transformations (5.62) and (5.63), this procedure generates a family 
{{lAm{y)^Vm{'y),y^m{y), Zm{y))}m& of rational solutions of the Painleve-II system (5.19). At the level of 



the logarithmic derivatives given by (5.38) we have a family of rational solutions of inhomogeneous Painleve- 



II equations with certain quantized values for the inhoniogeneity parameter a. In fact, it is known 22J that 
the equation 

V"{y) - 2V{y f + '^-yV{y) - (5.64) 

has a unique rational solution V — Va{y) if and only if a G Z. These may be represented as logarithmic 
derivatives of other rational functions Ua that are ratios of consecutive Yablonskii-Vorob'ev polynomials 
|13L 114) . which are known to have a number of remarkable properties. 

It is important for us to recognize that the poles of these solutions can occur for finite real values of y. 
For example, applying (5.62) twice with Wo(j/) = 1 we obtain Ui{y) = —y/6 and l^2{y) = {y^ + 6)/(362/), and 
the latter has a pole at y = 0. These poles correspond to values of y for which the discrete Schlesinger maps 
from the neighboring values of m fail to exist, and consequently for these values of y and m there exists no 



solution to Riemann-Hilbert Problem 4-1 



While there exist values of y for which the Backlund transformations (5.62) or (5.63) may not make sense. 



we may always interpret these transformations as maps on the ring of rational functions, and as such they 
can be shown to have a kind of singularity confinement property: if yo G C is a pole of the solution for some 
m, then yo is a regular point of the solution for mil. Indeed, a local analysis of the system (5.20) in the 
spirit of Painleve's original method [19 shows that all poles of {Um{y),^m{y)) are necessarily simple and 
simultaneously occur in both functions; if yo is a pole, then for some constants fc G C \ {0} and uj £ C we 
necessarily have 

"1 yo, . / 1 



^rn (y) = k 



y-yo 



18' 



+ ^{y-yo)+ iy~yor + Usiy) 



12 



V,„(y) = -fc" 



1 



ya 



(y-ya) 



1 

12 



{y-yof + V3iy) 



(5.65) 



.y-yo 18 

where Uz{y) and Vj,{y) are analytic functions each vanishing to third order at y = yo. Then, using these 
formulae in the Backlund transformations ( |5.62 )-(5.63) shows that the singularities at y = yo in the pairs 
{Um+iiy),Vm+i{y)) and {Um-iiy),Vm-i{y)) are aU removable. Moreover, both V.,„+i(y) a,ndUm-i{y) have 
simple zeros at y = yo- It is easy to see that the singularity confinenment property extends to the Schlesinger 
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transformations that generated (5.62)-(5.63). That is, if Zm{C',yo) does not exist (because yo is a pole of 
l^m{y) and Vm(y)), then both Zm+i(C; yo) and Z„j_i(C; j/o) do exist and they may be calculated from Zm(C; y) 
by applying Schlesinger transformations and then taking the limit y ^ yo- We summarize these results in 
the following proposition. 



Proposition 5.2. Let m £ Z. 



or (5.59) applied to the base case of m 



Then by iterated Schlesinger transformations given by (5.531 subject to (5.57) 
characterized by Proposition 5.L a matrix function ZimiC'jy) 
well-defined as a rational function of y with simple poles. The poles of Zm(<^; y) are exactly the poles of the 
corresponding potentials {Um{y),Vmiy)) obtained from the base case of {Uo{y),Vo{y)) = (1,— ?//6) by iterated 
Bdcklund transformations ( |5.62 |-(5.63). 

• If 2/0 "is not a pole, then there is a neighborhood A (Z C with yo £ A such that for y £ A, the matrix 

^m{C;y) 

— is analytic in y, 

— is the solution of Riemann-Hilbert Problem\4.1\ 



however both Zm+i{C] y) 
liCiy) fl'^c regular at yg (in the sense of taking limits from y ^ yo), and both Vm+i{y) o,nd 



yo; 



^ has an expansion for large Q of the form (5.1 1 differentiable with respect to both y and ( and 

uniform for y G A, and 
— is uniformly bounded for y £ A and C in compact subsets of C 
• Ifyo is a pole, then Riemann-Hilbert Problem 4-1 has no solution for y 
and Z 

lAm-i{y) have simple zeros at y — yo- 
Since the poles are simple, the bound on the elements o/Zm(^;?/) (and its inverse, since det(Z„i(C; — 1 
where defined) can be strengthened to 0{\y — 3^{hlrn)\~^) (or equivalently 0{\y — 3^{Vm)\^^)) uniformly on 
compact sets in the C,-plane. This improved bound also applies to all of the coefficients in the expansion 
(5.1) and to the error term as well. In particular, this implies that at y — yo £ 3^(Um) = ^{Vm), 
the functions Hm{y) = \Am,22{y) = -\Am,ii[y), Bm,i2{y) = 2H„iiy)U,niy) - U[^{y), and B,„,2i(i/) = 
^m{y) ~ '^Hm{yWm{y) all have simple poles. 



6. Inadequacy of the Global Model 



On one hand, it seems quite reasonable to expect that the global model Pm{w) defined by (4.66) should 



provide a good approximation to P(w), at least if the integer m is properly chosen. Indeed, we have shown in 
an earlier paper [7] that far from criticality a global model obtained by steps completely analogous to those 
we have followed here (choice of an appropriate g-function, opening of lenses, pointwise approximation of 
jump matrices away from points of nonuniformity and the solution of a corresponding outer model problem, 
and exact solution of the jump conditions near points of nonuniformity yielding inner models that match 
well onto the outer model at disk boundaries) leads directly to a Riemann-Hilbert problem for the error of 
small- norm type. On the other hand, in the present case the formally large exponents is/e^r were removed 
from the jump matrices in our construction of the inner model, and we may expect them to reappear when 
the mismatch between the inner and outer models is calculated on dU, contaminating the error estimates 
when s is small but large compared with ejy. 

In order to properly gauge our prospects for success, we introduce a new matrix unknown Qm(w) con- 
structed from the old unknown P{w) and the explicit global model Pm{w) defined by (4.661 as follows: 

Q„(u;) :=P(?«)P„(u.)-i. (6.1) 

We might expect to be able to prove that Qm(w) is a small perturbation of the identity matrix with the use 
of small-norm theory applied to the Riemann-Hilbert problem satisfied by this matrix function. 

To begin to determine the nature of the Riemann-Hilbert problem satisfied by Qm{w), firstly note that 
Qm{w) is analytic for w £ U, since for such w, P„i{w) = Pml""^) determinant one and satisfies exactly the 
same jump conditions as does PJw). For w £ C\U, Qmiw) is analytic except on the arcs of the jump contour 
for P{w) illustrated in Figure 



3.2 



it can further be checked that since both P{w) and P,„(w) = P°^*^{w) 
change sign across the contour segment (b,0), Qm(w) also extends analytically to this segment. The jump 
contour for Qm(iy) also contains dU, where P{w) is analytic but Pmiw) has a jump discontinuity stemming 
from the mismatch (because C is large but finite on dU) between P^)^{w) and P™*(it;). The jump contour 
for Qm(w) is illustrated in Figure 6.1 
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Figure 6.1. The contour E of discontinuity of the sectionally analytic function Qm(u'). 
With the exception of the circle dU which is clockwise- oriented as shown, all other arcs are 
oriented exactly as components of the jump contour ofF{w). 



Let S denote the jump contour for Qmiw), that is, Qm(w) is analytic for w G C \ S, and it takes its 
boundary values in the classical sense on this contour from each component of the domain of analyticity. 
Also, since both P(w;) and Prn{w) tend to the identity matrix as w ^ cx) (the former by hypothesis and the 
latter by the exphcit formula for P™*(u;)), we have Qm(w) — >■ I as w — >■ oo. 

For fixed m G Z, it turns out that for {x, t) sufficiently close to criticality, Q,n+(C) — Qm-(C)(I + ^{^n)) 
as Eat — )■ for ^ e S \ (1R+ U dU), with the estimate being uniform with respect to ^ on the specified contour 
arcs and {x,t) near criticality. Indeed, exactly at criticality we have q(w) = g{w) where g{w) = g{w;x,t) is 
as defined in [71 Section 4] and so by continuity with respect to x and t we will have near criticality that the 
jump matrix for 0{w) is subject to the same error estimates as recorded in [71 Proposition 5.1] generally 
valid away from criticality, where the neighborhoods Ui and U2 in the statement of the proposition are taken 
to coincide with U . Making the transformation from 0(u;) to P(?x;) we learn that the jump condition for 
P(w) is of the form P+(0 = P-(0(I + C'(ew)) uniformly for all ^ in the jump contour pictured in Figure 3.2 
with the exception of ^ € J7 and on the ray ^ > —1. But, Pm(w) satisfies exactly the same jump conditions 
as does 'P(w) within U and in the interval — 1 < ^ < 0. The estimate on the jump conditions for Qm(w) 
asserted at the beginning of this paragraph then follow by exact computation. 

For ^ S K+ \ / (recall that / is a small open interval containing ^ = 1 as shown in Figure 1.6 1, both P(w) 
and Vrniw) = P°"*(w) satisfy the same jump condition, from which it follows that for ^ > but ^ ^ /, 

Q„,+ (0 = V+{i)VZX{i)-' = cT-2P-{S)a2 [a2P°r (0^2] = a2P_(e)Pr-(0"V2 = a2Q™-(e)a2. (6.2) 

For ^ G /, it follows from the fact that 0+(^) +fl-(0 = for ^ > that |7i Proposition 3.3] holds, implying 
that P(it;) satisfies P+(C) = tT2P-(0<''2(I + C(e7v))- On the other hand, by construction Pm{w) satisfies 
the corresponding jump condition without error term. Because Pm(w) is uniformly bounded with bounded 
inverse, by a calculation similar to (6.2| we see that Qm+(0 = o'2Qto- (^)cr2(I + C'(eAf)) holds uniformly for 

C e / c M+. 

To compute the jump condition for Qm(w) across the negatively-oriented boundary dU of the disk U , 
first use the fact that P(w) is continuous across dU to obtain 



Q™+(o = p(Opr(o~' = p(c)pi^(e)-' pj^(Opr(o~' =q™-(o vtmz'io-' 



By substitution from (4.601 and (4.65), we see that 



_^(2m-l)(T3/6 iso-3/ejv 



??rn(e)"^'Z,„(C;y)(-C)(^-'"^^^/'r7™(e)"'^^l e-5-3/^«,a-2™).3/6^^^^) 



e e dU. 
(6.3) 



(6.4) 
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for ^ G dU, where C = e^^^WiO- 

We therefore see that for {x,t) close enough to criticahty, Qm{w) satisfies the conditions of the following 
type of Riemann-Hilbert problem. 

Riemann-Hilbert Problem 6.1. Seek a matrix Qm(w) with the following properties: 

Analyticity: Q„i(if) is analytic for w e C\S, where S is the contour (independent ofx, t, m, and 



^n) pictured in Figure 6.1 and in each component of the domain of analyticity is uniformly Holder 
continuous up to the boundary. 

Jump conditions: The boundary values taken by Qm(u') on S are related as follows. For ^ G 
dU C S with clockwise orientation, 

Qrn+(e) = Qr«-(e)VQ,„(C) (6.5) 

with jump matrix Vq^^ (^) defined by the right-hand side of ( |6.4[ ). For ^ G M+, 

Qm+(0 - 'T2Q„-(0'T2(I + X„,(0) (6.6) 

where Xm(^) = ©(e^v) holds uniformly for ^ G / and Xm(^) — for ^ G M_|- \ /. For all remaining 
^ G S we have the uniform estimate 

Q„+(0 = Qm-(0(I + (^(ejv)) (6.7) 

as ejv i 0. 

Normalization: The matrix Qm(w) satisfies the condition 

lim Q„(w)=I. (6.8) 

I arg(-t«)|<7r 

of the above jump matrix estimates are uniform for bounded m. 
Note that since the jump matrix for (^m{w) is not specified here precisely, the conditions of Riemann- 



Hilbert Problem 6.1 are not sufficient to determine Qm(t«), but we shall see that they are enough to approx- 
imate Qm(w) with suitable accuracy as e^v ~^ 0. 

Now VF(^) is bounded away from zero for ^ G dU , so from the normalization condition on Zm(C;y) given 



as part of Riemann-Hilbert Problem 4.1 and understood in the stronger sense of (5.11, and from the fact 



that ??m(0 its reciprocal are uniformly bounded as e^r — >■ for ^ G dU , we obtain 



Pm(OPr(0"' = ei^™"'^"'^^e5^'^^/^" l + 0{e'j^') e-5*"^/^"e^"""'^"^'°dw(e)-"% ^€dU. (6.9) 



1/3 



<^-M- (6.10) 



Recall that d]\[{w) = 1 + ©(e^r) uniformly for w G [/. It follows that the jump matrix for Qmiw) will be a 
small perturbation (of order O(e^) for some ii G (0, 1/3)) of the identity for ^ G dU if 

2m — 1 s 
3 eArlog(e^^) 

Subject to this condition, we may formulate a Riemann-Hilbert problem for Qmiw) of small-norm type, that 
is, a problem that can be solved by iteration. This calculation leads to the suggestion that different choices 
of the integer m will be required to control the error in different regions of the (y, s) parameter plane, an 
idea that has proven useful in at least two very recent works on related problems of semiclassical analysis 
for nonlinear waves. In the work of Claeys and Grava [11^ on the small-dispersion limit of the KdV equation 
near the "trailing edge" of the oscillatory wavepacket generated by a gradient catastrophe in the solution of 
the limiting Burgers equation, the authors find that different choices of an integer parametrizing an outer 
model problem are required as the rescaled distance from the trailing edge varies. Similarly, in the work of 
Bertola and Tovbis |4] the semiclassical solution of the focusing NLS equation is examined for points near 
the "primary caustic" or "breaking curve" generated from an elliptic umbilic catastrophe in the solution of 
the formal limiting Whitham system, and it is again found that different choices of parametrices must be 
made as a certain rescaled distance from the caustic curve varies. In both of these problems, as will be the 
case here, the rescaling of distance involves logarithms of the small parameter to get the balance right, and 
the resulting approximate solutions have a "solitonic" nature (in the KdV problem they are simple solitons 
while in the NLS problem they are soliton superpositions forming breathers). 
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To implement this suggestion, we will choose a value for m and then define corresponding regions of the 
{y, s) plane in which we can make use of this value of m. Also, we will find it necessary to relax the small- 
norm condition (6.101 in order to cover a large enough part of the (j/, s) plane with overlapping regions, and 
this means that the Riemann-Hilbert problem solved by the matrix Q,„(ti;) will no longer be treatable by 
iteration, and instead we shall have to improve our inadequate model Pm(w) by building parametrices for 
Qm{w), another technique used in the papers [TT] and [3] (also in the more recent work [S] of Bertola and 
Tovbis where the solution is analyzed in a neighborhood of the elliptic umbilic catastrophe point itself) . 

For each integer m G Z we define a shifted and scaled coordinate p,„ equivalent to s by the relation 



2™ . / 
s = — eArlog(e^ ) + ejvp„ 



(6.11) 



Definition 6.1 (Sets il^). Let m For suitably small 6 > and arbitrary k> 0, set 



and 



(y, s) G with \pm\ ^ - logC^Tv""^) such that 
3 

\y- ^{U,n)\>6e--P- 
\y-^illrn)\<Se-iP- 



{y, s) £ with < - log(ejY"'^) + k such that 

o 



«/ - 2 log(eAr^) - ii<Prn <0 and 
if < p,n < ^ log(e^^) + k| 



(6.12) 



\y - ^(V™)| > ,5e-5P— 1 1/ < pm-i < 3 log(e^') + k and 



\y - ^iVm)\ < 6e^P"^-' if 



log(e^^) - K < Pm-i < 



(6.13) 



The set 17+ is basically a horizontal strip with excluded notches located near the poles of Urn{y) and 
upward-extending teeth located near the zeros oiUmiy)- Similarly, fi" is a horizontal strip with excluded 
notches located near the poles of Vm{y) and downward-extending teeth located near the zeros of Vm(y). 
Since the poles of Um and Vm coincide (although their zeros do not in general), that is, S^iUm) = ^(Vm), 
the notches of 11+ align with those of 17" . This is very important, since we will make use of the solution 
Zm(C;y) of Riemann-Hilbert Problem 4.1 for {y,s) £ 17+ U 17^, and therefore we must exclude from this 



union points with y-values too close to points of £^(JJ„i) ~ £^{Vm) at which ZmiC'^y) fails to exist. See 
Figure |6.2[ 

Since Vm+iiy) — 1/Wm(2/) according to (5.621, it is clear that i^{Um+i) = SfiVm) and ^{Um+i) = 



^O^m)- It follows that while the sets 17+ and i^„^^i have disjoint interiors, their union is a horizontal strip: 



17+ U 17, 



m+l 



{y, s) e such that \pm\ < 7^ ^os{^n) + « 



(6.14) 



The two regions interlock near the singularities of log \Um (y) \ (equivalently the singularities of log | Vm+i iy)\)', 
the teeth of 17+ fit in a one-to-one fashion into the notches of 17,~^ ]^, and the teeth of 17^^^^ i^- a one-to- 
one fashion into the notches of 17+ . See Figure 6.3 In Figure 6.4 we illustrate the tiling of a region of the 
(y, s)-plane with many such strips when k = (the strips overlap if k > 0). This latter figure illustrates 
another feature of the k = tiling: since by the Backlund transformations ( |5.62 |-(5.63l, yo is a simple zero 



of Ufn-iiy) if and only if j/g is a simple pole of both Umiy) and Vmiy) if and only if yg is a simple zero of 
Vm+i{y), the tips of the teeth of ^^-1 matched in a one-to-one fashion onto the tips of the teeth of 

17; 



m-t-l ■ 



We will use Definition 6.1 in the following context: we fix a bound B > and consider only indices 
TO €E Z with |to| < B; we then choose S > sufficiently small that for all to e Z with |to| < B, no two real 
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zeros of Um{y) 



F 




, ^ . / ^ . / , 






' 

zeros of Vm{y) 








poles of Um{y), equivalently, poles of Vm{y) 



Figure 6.2. The regions f2^j (pink) and f2,~ (blue) for m = 6. For this figure, k = 0, 
S — 0.2 and ejv — 0.1. Note the alignment of the notches. It is for (y, s) in the union of 
these regions that we select the value m — 6 in both the solution P^*'(w) of the outer model 
problem and the solution Z„(C;?/) of the inner model problem. 



S 



i 


L 

poles of Vm+i{y), equivalently, zeros of Um.{y) 


/ \ / \ 1 \ 




poles oiUrniy), equivalently, zeros ofVm+i{y) 




y 



Figure 6.3. The regions 51+ (pink) and r2„j^]^ (blue) for m = 6. Their union is a closed 
horizontal strip in the (y, s)-plane. For this figure, k ^ 0, S = 0.2 and ejv = 0.1. 



singularities of log \Urn{y)\ or log |Vm(y)| are closeij^ than 26. This prevents neighboring teeth/notches from 
overlapping at their widest point (corresponding to = 0). 

6.1. Parametrix for Qm(w) in ^ti- Even with the choice of the integer m corresponding to {y, s) € £7^ it 
turns out that Riemann-Hilbert Problem |6.1| cannot generally be solved with the use of small-norm theory 
as the jump matrix Vq^(^) is not close to the identity for ^ e dU, at least not uniformly for (y, s) G 51^. 

^Preliminary asymptotic analysis suggests that 5 = 0{B~^^^) as iJ i~ oo, which is one indication of failure of uniformity 
of the analysis in this paper when t ^ s 2> e^f log{eJ^^). Indeed, for t > fixed and x x^rit fixed one expects a different 
description of the sine-Gordon dynamics for sufficiently small ejv in terms of higher-genus Riemann theta functions. 
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^eivlog(e^^) 



14 











1 1 1 1 1 



-10 







10 



Figure 6.4. The regions riig, . . . , riy (pink) and ri_g, . . . , (blue), for the "tiling" case 
K = 0. 



Therefore, even though we may think of Clm{w) as an "error" matrix gauging the mismatch between the 
matrix function Fjn{w) and the matrix function P{w) for which the former was introduced as a model, we 
must in general carry out further modeling of Qm(w) by creating an explicit model (parametrix) for it, and 
comparing the model with Qm(w) iu hopes of — finally — arriving at a small- norm problem for the matrix 
ratio. We now proceed to build such a parametrix for the error by a procedure referred to in some papers 
|10|, lllj as improvement of the parametrices. 

Let TO G Z with |to| < i? be fixed, and suppose that (y, s) G $1+ with y G M bounded. In particular, 
this implies that the coordinate p,„ defined by (6.11 ) is subjected to the inequalities \pm\ < | log(£jv^) + 
Using the coordinate Pm, and recalling Proposition 5.2 along with the expansion (5.1), we write the matrix 



Vq^(^) defined by the right-hand side of (6.4 1 in the form (assuming only that y is bounded) 



Vq„(C) = exp ( ^ (p,n + ^ log(e/) ) erg 



exp ^-^ + ^ log(e^i)^ 0-3^ dN{0~ 



VmiO-''' (6.15) 



because W{^) is bounded away from zero on dU. Since r/„i{^), (iAr(^), and their reciprocals are all bounded 
on dU and Pm > ~l log(e^^) - k for (y, s) G 17+ , 



Vq,„(0 = exp ( ^ (pm + ^ log(e/) ) erg 



exp 



2/3 n„ 



|y-^(Z^m)| 



(6.16) 



holds for {y,s) G fi+ with y bounded and ^ G Now, we want to neglect all but the (l,2)-entry of 

the matrix e](^ Am{y)W{^)~^ + e^'^Bm(j/)M^(^)^^. Since the matrix elements of Am{y) and Bm(y) are 

0{\y— I^{Um)\~^), we may do this at the cost of adding a new error term of the form 0{e]^'^\y— ^^(Um)\^^) 
(dominated by the diagonal entries). This error term dominates the one already present in ( |6.16 1 for 
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{y, s) € r2+ as a consequence of the inequality p„i < | log(e^"'^) + n. Therefore, 

/ 1/3 \ 



Vq„(0 
where, defining 



1 Y+iC) 
1 



-TV 



the expUcit function appearing in (|6.17|) is 



e e 9C/, (y, s) £ n+ with y bounded, (6.17) 

(6.18) 

(6.19) 



H := eP'-h+iw) (^A,r,MyW{w)-' + e](^ B„,^,2{y)W{wy 
= eP-h+{w) (Urn{y)W{wr'+e]^'[2H^{y)U^{y) ~U:MW{wy 

Note that Yj^{w) extends from ^ e dU as a meromorphic function of w € U whose only singularity is a 
double pole at w = and that Y^{(^) also depends parametrically on pm and y. Since div(0 = 1 + C(eAf) 
holds uniformly for ^ G another way to write (6.17) is 

-1 / 



1 YrtiO 

1 



= 1 + 



1/3 

Iy-^(Z^™)I 



,1/3 



\y-^{U,n)\J' 

C e af7, {y, s) e n+ with y bounded. (6.20) 



Proposition 6.1. The following three estimates hold uniformly for ^ G dU and y bounded. 
In the upward-pointing "teeth" of ^V^: 



1 y;ii{0 

1 



= I + 0(e]v^'eP'"|y-^(Z^„)|), 

(y,s)er!+, p,n>Qand\y- ^{U^)\>5e-P-^. (6.21) 
/n f/ie pari of resembling a "strip with notches": 



Vq,„(0 
Elsewhere in f2+ 



1 y,t{0 

1 



,1/3 

-AT 



|y-^(Z^,n)l 



1 Y+{0 
1 



(6.22) 



(6.23) 



Proof. First suppose that < p„i < |log(£Ar^) + (so we are considering (y, s) to lie in the upward- 
pointing "teeth" of 51+ ) and that y is bounded. In this sub-region, \y — 3^{Um)\ is bounded away from zero 
because the teeth are localized near the zeros 2f{Um) oiUm{y). Therefore 2HjnUm ~ ^'n is bounded and 
Um[y) = 0{\y - ^ (^m)|), so Y+(C) = 0{eP"-\y - ^{Uyn)\) + 0(l) holds uniformly for ^ e dU . Putting this 
estimate into ( |6.20[ ) gives 



1 



I + 0{^^'') + 0(e]v/'eP-|y - ir(Z^„0|), 



i e 9?7, (y, s) e 51+ with y bounded and p„i > 0. (6.24) 



> 0. 



This proves the estimate (6.21 ), and it also proves the estimate (6.23) in the case that p„ 

'|log(e^^) — K < pm < (so we are considering (y, s) to lie in the com- 
■'strip with notches"), and that y is bounded. Since according to 
i{y) —U'^{y) have simple poles at the points of ^{Um), it follows 
- ^{Um)\^^)- But in this subre gion we also have the inequality 
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Now suppose instead that 
plementary part of f2+ that resembles a 



Proposition 5.2 hothUmiy) ^^'^ '2Hm{y)Ur, 
that in this sub-region, y,^(0 ~ 0{eJ'"^\y - 



\y — 3^{Um)\ > (5e2P™, so with p,„ < in fact we have Y^{(_) = 0(1). Putting this information into (6.20) 
gives 



Vq(0 



1 Y+iO 
1 



o 



^1/3 



, ^ £ dU, {y, s) e fi+ with y bounded and p,n < 0. 



This proves (6.22), and it also proves (6.231 in the case that p„i < 0. 

Corollary 6.1. The following estimate holds uniformly for ^ G dU and {y, s) G f2^j with y bounded: 



VQ„(e) 



1 y^tio 

1 



(6.25) 

□ 



(6.26) 



While uniform, this estimate fails to be sharp except possibly at the vertical extremes of where \pm \ ~ 
^log(e^'). 

As Eat J, 0, the function Yj^{^) is not uniformly small on dU as (y, s) ranges over £7+ with y bounded. 
This implies that to approximate Qm{w) for {y,s) e it will be necessary to build a parametrix for it 
that takes into account the explicit terms in (6.17). 

Riemann-Hilbert Problem 6.2 (Parametrix for Qm(w) in ^ti)- ^^'^^ numbers y andpm be given, and 
assume that y ^ .^(Um)- Seek a matrix Q^(w) with the following properties: 

Analyticity: (w) is analytic for w G C \ (M+ U dU), and is uniformly Holder continuous up to 
the boundary of its domain of analyticity. 

Jump conditions: The boundary values taken by (w) on the contour U dU are related as 
follows: 

Q™+(C)=^2Q+-(0^2, C>0, (6.27) 

and 

[1 Y+iO] 



1 



where dU is taken to be negatively (clockwise) oriented and (0 is defined by (6.19). 
Normalization: the matrix Q,^(w) satisfies the condition 

lim Q+ («;)=!. 



(6.28) 



(6.29) 



We solve this problem uniquely for Q+ (w) as follows. Firstly, define a related unknown F+ (w) as follows: 

'Q+(w), w^U 



1 Y+{w) 
1 



w €U. 



(6.30) 



Clearly, F+ (w) is continuous across the contour dU. Now, rjm{w) is analytic in [/, but W{w) has a simple 
zero at w = ~ —1, so F+ (w) generally has a double pole ai w = w■^,. To handle the only remaining jump 
condition for F+(w), namely F^j_|_(^) = o'2F^j_ (^)cr2 for ^ > 0, we introduce r := z(— w)^/^ and set 



G+(t) : = 



F+(r2), 3{r} > 

a2F+(r2)a2, 3{r}<0. 



(6.31) 



It is now evident that G^j(t) is a rational function of r G C with simple poles at r = ±iS where S := 
(— w*)^/^ > 0, normalized to the identity at r = oo, and satisfying the symmetry relation G+(— r) = 
cr2G+(r)cr2. 

To fully characterize the singularities of G+ (r), we will relate its principal part at r = iS to its regular 
part. Note that in a neighborhood of r = iS, Y^{t'^) has the Laurent expansion 



Y+{r') = , + + 0{1) 

[t — iby T — lb 



(6.32) 
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where 0(1) represents a function analytic at t — iS and where 

(6.33) 

Note also that J+ and are both real when y e M and p„j e M. Let us write the Laurent expansion of 
G^j(t) about T = iS* in the form 



G+(t) 



iS 



+ (e+ , f+) + isLK)ir - ^S) + 0{{t - iSf), 



(6.34) 



where 0((r — i'S)^) represents a matrix function analytic and vanishing to second order at r = iS. Since 
Q+(w) is meant to be analytic at w = it;*, from (6.30) we see that the principal part of G+(r) at r = i5 is 
determined from its holomorphic part as follows. The first column vanishes: 

a+ = c+ - 0, (6.35) 

and the second column satisfies the linear relations 

b™ = ^me+ and d+ ^ J+g+ + iK+e+. (6.36) 

Taking into account the symmetry condition G^(— t) — cr2G^(T)cr2 and the notation of the principal part of 
the Laurent expansion (6.34), the matrix G^(r) necessarily can be expressed in partial fractions as follows: 

, (0,b+) , (0,d+) , (*a2b+,0) , (-*a2d+,0) 



G+(t) 



(6.37) 



(r-iS')2 T-iS {T + iS)^ T + iS ' 

where we have also used ( 6.35[ ). Evaluating the first column of this expression and its derivative at the 
regular point (for the first column) t = iS we find that, in the notation of the regular part of the Laurent 
expansion (6.34 ), 



4w* 



- and g+ = - 



25* 



^1 



1(72 



d+, s -.^ i-w.y/' > 0. 



(6.38) 



Substituting these expressions into the relations (6.36) yields a closed inhomogeneous linear system for the 
vectors b+ and di; : 



This system has a unique solution: 

4u;* J+ 



-cr2b+ 



i r 

47«, 



(72 b 



/+ 

2S 



(72 d+ = 



4w* 



(72 d 



7+ 


iK+ 





(6.39) 



b,1 
d.1 



\w,{l%wl + 2SJ+K++i[J+Y) 
D+ [ [J+f -lQwlJ+-Z2wlSK+ 

4*5 r 4u;.(2[J+]3 + S[J+?K+ - l&wlSK+) 
D+ [[J+f + lQwl[J+f + i2wlSJ+K+ - mwl[K+f 

where the denominator 

:= {Wwl + [J+ff + Mwl{J+ + SK+f > 256wt 

is strictly positive, because w^, « —1. Substituting (6.40)-(6.41 ) into the partial fraction representation 
(6.37) completes the construction of G+(r). By writing r — i(— w)^/^ and restricting G+(r) to the upper 
half-plane we then recover F+ (w), and then from (6.30) we recover the solution Q+ (w) of Riemann-Hilbert 
Problem O 
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(6.40) 



(6.41) 



Proposition 6.2. For each fixed integer m with \m\ < B, Riemann-Hilbert Problem 6.2 has a unique and 
explicit solution (w) satisfying det(Q^(z«)) = 1 where defined. Also, Q^(w) is uniformly bounded for 
w d C \ (U U M+) and (jj, s) e 17+ with y bounded as e^r — > 0. 



Proof. The above construction clearly yields a solution of Riemann-Hilbert Problem |6.2[ and uniqueness 
follows easily from the conditions of that problem via a Liouville argument, which also shows that Q^(w) 
is unimodular. Since Q+(u') = F^(w) for w € C \ U , the uniform bound on Q+('«;) will follow from a 
corresponding bound on F^(w) valid when w is bounded away from the singularity at w*, or what is the 
same, a bound on G^(r) valid when r is bounded away from both iS and —iS where S := {—w^,y^^ > 0. 



According to the partial-fraction representation (6.37) of (r), it is therefore enough to establish that the 
elements of the vectors b+ and d+ remain bounded as e^r J, uniformly with respect to (y, s) G il^ with 
y bounded. But from (6.40) we see that the components of b+ and d+ really only depend on J+ and isr+ 
(since ~ —1) and while the latter can become unbounded, it is still obvious that b+ and d+ cannot. □ 

6.2. Parametrix for Clm{w) in il^j. Fix an integer m with |m| < B, and suppose that (y, s) G fl^ with 
y G M bounded. By introducing the scaled and shifted coordinate Pm-i defined by (6.11) and analyzing 



the jump matrix Vq^^(^) for ^ G dU as in the beginning of 5 6.1 we see that in the current situation a 

(6.42) 



lower-triangular model for the jump matrix is more appropriate. In particular, defining 
and then 

Y-,{w) := e-P-^h-{w) (v.m{y)W{wr' + e]i^ [V^v) ^ 2H„,{y)V,n{y)] W{w)-^) , (6.43) 

(note that (^) extends from ^ G dU to a meromorphic function for ^ G C/ whose only singularity is a 
double pole at w = w^,) we arrive at the following. 

Proposition 6.3. The following three estimates hold uniformly for ^ G dU and y bounded. 
In the downward-pointing "teeth" offl^: 



Vq„(0 



1 



= l + 0{e]i'e-P-^y~ 2f{V,n)\), 



{y,s) G n„„ 

In the part of fl^-^ resembling a "strip with notches": 

,1/3 

.(0 " =i + o' 



Vr 



1 



\y-^iVm)\ 



p„,_i<0and\y-3r{V,„)\>5eP" 



(y,s)Gf7-, \y-^iV„,)\<5. 



Elsewhere in fl„ 



Vq„(0 



1 

Y-{0 



Oie 



l/3x 
N ) 



(6.44) 



(6.45) 



(6.46) 



The proof of this statement is virtually the same as that of Proposition 6.1 The analogue for i7„ of 
Corollary |6.1 1 is then this result. 



(6.47) 



Corollary 6.2. The following estimate holds uniformly for ^ G dU and (y, s) G r2„j with y bounded. 

-1 



Vq„(0 



1 

While uniform, this estimate fails to be sharp except possibly at the vertical extremes of fl^ where |p,„_i| ~ 
|log(e^'). 

Based on this approximation result for Vq^(^), we build a parametrix for Qm(w) to satisfy the following 
problem. 

Riemann-Hilbert Problem 6.3 (Parametrix for Qm(w) in 12^). Let real numbers y and p„i-i be given, 
and assume that y ^ ^{Vm)- Seek a matrix Q^(w) with the following properties: 
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Analyticity: Q~ (w) is analytic for w G C \ (M+ U dU), and is uniformly Holder continuous up to 
the boundary of its domain of analyticity. 

Jump conditions: The boundary values taken by Q^iw) on the contour M_|_ U dU are related as 
follows: 

Q™+(0=^2Q--(0^2, ^>0, (6.48) 

and 

1 0" 



(6.49) 



where dU is taken to be negatively (clockwise) oriented and (^) is defined by (6.43) 
Normalization: the matrix Q~ (w) satisfies the condition 

lim Q-(u;)=I. 



(6.50) 



This problem is solved completely analogously to Riemann-Hilbert Problem |6.2[ Namely, we introduce a 
matrix G^{t) by 

'f,- (t2), Q{t} > 

^asF,- (r2)a2, 5{t} < 0, 

where we define F^{w) by 

■Q„(w), w^U 



Gr„(T) 



(6.51) 



Q,„(w) 



1 

Y-Aw) 1 



w eU. 



(6.52) 



Then, G~ (r) is a rational function with G^^{oo) — I and double poles at r = ±iS where S := (— w*)^/^ > 
only, satisfying the symmetry relation G^j(— r) = (T2G^j(T)cr2. 

The singularities of G~(t) are characterized by first noting the Laurent expansion of Yj^{t'^) about 
T = iS: 



[T-isy 



iS 



-0(l) 



(6.53) 



where 0(1) denotes a function analytic at t = iS* and where 

J . . . ^ , , I 111 . \C ^ 



e~^--^h'^{u,^^.^ (2i/„(y)V,„(y)-V;„(y) 



2w,h-^{w,)W'{w,fV^{y) 



4w,SW'{w,)^ 

+ e]/'(2H,n{y)Vrn{y) - vUy)) ( - 2w,h;:{w,)w'{w,) + 2w,h;MW"{w,) + h;MW'{w, 

(6.54) 

are real- valued coefficients when both y € M and Pm-i G The Laurent series of G~ (r) about t = iS has 
the form 

Gn^ir) - ^p^^ + ^^^^ + (e,- , f-J + (g,- , h-J(r - ^S) + 0{{t - zSf), (6.55) 
(r — i5)^ T ~ lb 

where 0{{t — iS)"^) represents an analytic function vanishing to second order at t — iS. Since Qm(u') is 
meant to be analytic at w = if*, from (6.52 ) we learn that the principal part of the second column of G^(t) 
vanishes: 

(6.56) 
(6.57) 



Ki = d,„ = 0, 



and the principal part of the first column satisfies the relations 
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Using the symmetry condition G„j(— r) = cr2G„(T)cr2 and taking into account (6.56 1, we see that the partial 
fraction expansion of G^{t) necessarily has the form: 

, (am,0) (c„,0) (0,-icr2a„) , (0,icr2C-) 



G™(r) 



iS 



(T + isy 



iS 



(6.58) 



By evaluating the second column of this expression and its derivative at r = iS* we can express and h„ 
in terms of aZ, and cZ, : 



1(72 



^2 

2S' 



and h„ 



0-2 



1(72 



5:= >0 



(6.59) 



Combining this with (6.57) yields a closed inhomogeneous linear system of equations for the components of 
the vectors a~ and c~: 



As before, this system has a unique solution: 



-0-2 









(6.60) 



4w* J,„ 



4i5 



16^2 j„ - 32w25if^ 
4w;,(16«;2 + 25J-if- + 3[J-]2)J 

[J-]4 + 16u;2[J^]2 + i2wlSJ-K- - i2wl[K- 
4u;,(2[J-]3 + S[J~YK- - l&wlSK-) 



where the denominator 



-.^ [IQwl + [J-J2)2 + 64u;2( J,- + > 2hQwt 



(6.61) 



(6.62) 



— 1. This essentially completes the construction of Q,„(w) solving 



The analogue of Proposition 6.2 in this case is the following. (Its proof is also nearly the same.) 



is again strictly positive because w» 
Riemann-Hilbert Problem 16.31 



Proposition 6.4. For each fixed integer m with \m\ < B, Riemann-Hilbert Problem 6.3 has a unique and 
explicit solution Q~(w) satisfying det(Q~(w))) = 1 where defined. Also, Q~(w) is uniformly bounded for 
w G C \ (f/ U M+) and {y, s) G il^j with y bounded as en — > 0. 

6.3. Accuracy of the parametrices for Qm{w). Error analysis. The error in approximating Qm{w) 
by the parametrix Q^j(w) is the matrix E^(u;) := Qm(w)Q^(w)^^. Since the jump contour dU for Cl^^{w) 
is a subset of the jump contour E of Qm(w), and since det(Q^(i(j)) = 1 for w G C \ dU, it is clear that 
E^(k;) is, like Qm(w), analytic for w G C \ E, and takes its boundary values on E in the same classical sense 
as we require of Qm{w). Also, E^j(w) — > I as w — ?> oo since this is true of both factors. 

We now consider the jump conditions satisfied by E^{w) as a consequence of those known to be satisfied 
by Qm(w) from the statement of Riemann-Hilbert Problem 6.1 and those satisfied by Q^(u') from either 
Riemann-Hilbert Problem |6.2| or |6.3[ If ^ G , then 



= ^TaQ™- (0^2(1 + X„,(e))a2Q±_(0" ^2 



(6.63) 



a2E^_(e)a2 • <72Qi_{0<72iI + X,„(0)a2Q,=f,_(0"V2 



If we write the matrix in brackets above as I -I- X^(^), we see that X^(^) = wherever Xm(^) = (that 
is, everywhere except in the interval / near ^ = 1). Using Proposition 6.2 and the fact that Xm(^) = ©(ejv) 
for ^ G /, we see that under the condition that {y,s) G il^j with y bounded we also have X^(^) = ©(cat) for 



^ G /. Similarly, Proposition 6.4 shows that under the condition that (y, s) G r2„j with y bounded the same 

(6.64) 



estimate holds for X^ (^) . If ^ G E \ (M+ U dU) , then the parametrix has no discontinuity, so 

E^.+(e) = Q™+(OQ^.(0"' =E±_(OQ^„(0(i + o(eA^))Q?„(0"' 
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Therefore according to Proposition 
C'(eAr)) for such ^. Analogously, 



6.2 



if (y, s) _e_f2+ with y bounded, we will have E+ + (^) = E+_(^)(I + 
shows that {y, s) £ f2~ with y bounded implies that 



'roposition 



6.4 



E- , (0 = E;^_(e)(I + ©(ejv)) for such ^ Finally, if ^ G 5f/, then 



E++(e) = Q™+(OQ++(C)" 
= Q™-(OVq,„(0 



1 



E+-(e) Q+-(?)Vq„(0 



1 



(6.65) 



If (y, s) G rj^j with y bounded, then according to Propositions 



6.1 



and 



6.2 



we find E+(0=E+_ (0(1+0(1)) 



where the size of the error term is 0{e]^^) at the worst (according to Corollary 6.1) but is more typically 
smaller according to the more precise estimates (6.21 1-( [6.23 1 enumerated in Proposition 6.1 Similarly, for 

-1 



E- + (e) = Qr: 



1 



E_(0 Q-_(C)Vq„(C) 




(6.66) 



E„_(e)(i + o(i)) 



so if (y, s) G fi^ with y bounded, then from Propositions 
with a similar characterization of the error term. 

The dominant terms in the deviation of the jump matrix from the identity come from ^ e dU. Therefore, 
we observe that the matrix E^ satisfies the conditions of the following type of Riemann-Hilbert problem. 

Riemann-Hilbert Problem 6.4 (Small-Norm Problem for the Error). Let m be an integer with \m\ < B, 
and suppose that (y, s) £ il^j with y hounded. Seek a matrix E^(w) with the following properties: 

Analyticity: E^(u') is analytic for w G C \ S, and is uniformly Holder continuous up to the 

boundary of its domain of analyticity. 

Jump conditions: The boundary values taken by E^(ti;) on the contour S are related as follows: 

E± (0-a2E±_(0<T2VE±(0, e>0, (6.67) 



E± (0-e±_(C)Ve±(0, 



E^vs;, .--\-+, (6.68) 
where the following estimates hold for the jump matrix. Firstly, we have Vg± (^ ~ I for f > with 

I is compactly supported and vanishes identically in a neighborhood 

L-°(T,) = 0{e]l^), or, more precisely, ||Vg± - I||l~(s) = 



log(l^l) sufficiently large, so Vj 



o/ ^ — 0. Next, we have 
0{e^{y, s; Eat)), where for {y, s) G n+. 



e+(?/,s;eAr) 
and for (y, s) G 



^1/3 
Cat J 



p,r^>0and\y-3f{U^)\>Se-P- 
\y~^{U„^)\<6 



y-^{Vm)\, 



,iy,s;eN) := < e]y/'|y - ^(V™)| 



N 
^1/3 



elsewhere in f2^, 



p„,-i < and \y ~ 2f{V,n)\ > Se^" 

\y-^{V,n)\<S 

elsewhere in f2z, . 



Normalization: the matrix E^{w) satisfies the condition 

lim E±(«;)=I. 



(6.69) 



(6.70) 



(6.71) 
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Proposition 6.5. Suppose that m is an integer with \m\ < B , and that {y, s) G fJ^^ with y bounded. Then 
for sufficiently small ejv, there exists a unique solution E^{w) of Riemann-Hilbert Problem 6.4 The solution 
has expansions for small and large w of the form 




nE±](-Z«)-l/2 + 0(«;-l), «;^(X), 

and we have the estimates 

[°E±]=I + 0(e±(y,s;6^)), [?E±] = ^(e^ e^)), ^nd HE^J = 0(e± (y, s; (6.73) 



holding uniformly for {y,s) G with y bounded, where e^{y, s; epf) and e„(y,s;ejv) are defined by (6.69) 



and (6.70) respectively. 

Proof By the substitution F± (r) E± (r^) for 3{r} > and F± (t) := craE* (r2)cr2 for 3{t} < 0, we 



translate the conditions of Riemann-Hilbert Problem 6.4 into a list of conditions satisfied by F^(t). This 
"unfolding" of w to the r-plane implies that F^(t) is analytic on the complement of a compact contour S', 
along which it satisfies jump conditions of the form F^^^ = F,^j_Vp± where Vp± — I satisfies the uniform 
estimate ||Vp± — = 0{e^{y, s; en)) and e^{y, s; e^) is defined by (6.69) for {y,s) G f7+ and by 



(6.70) for (j/, s) G r2„. Since these estimates are all dominated by 0{e]^^), for sufficiently small e^r, the 



small- norm theory of matrix Riemann-Hilbert problems in, say, the L^-sense implies that F^j(r) may be 
constructed by iteration applied to a suitable singular integral equation, and it follows from that theory 
that an estimate of the form ||F^j — I||l°°(A') ~ ^{^miv^ ^n)) actually holds whenever K is a compact set 
disjoint from the jump contour S'. In particular, F^(r) is analytic on some compact neighborhood K of 
the origin t = 0, so it has a Taylor expansion that we write in the form 

Fi{T) = [°Fi] + [°Fi]r + Oir^), r 0, (6.74) 

and by the uniform estimate on F^(t) - I we have [qE.^J [qF^J = I + OdlF^^ - 1\\l=°(k))- Applying 
the Cauchy integral formula to F^(t) — I on a small loop in K surrounding the origin then shows that 
[jE^] = i[iF^] = OdjF^ — I||Loo(;f)). Also, F^(t) tends to the identity as r -> oo and it is analytic for 
sufficiently large |r| so it has a Laurent expansion of the form 

F^(T)=I+nF±]T-i+0(r-2), r^oo, (6.75) 



Now we can take K as an annulus enclosing a large circle. Applying the Cauchy integral formula on such a 

\l=°(k)) 



circle to F± (r) - I we find that ^E^^^J = -^HF.^J = 0{\\F^, - I||l-(k))- This completes the proof. □ 



The magnitude of the expansion coefficients of E+ (w) « I and of F^^_^_-f^{w) « I in the interlocking regions 



and respectively is illustrated in Figure 6.5 In particular it is clear that the order of accuracy 



1/3 

is 0{e^ ) except very close to the common boundaries between teeth and notches. Moreover, the order of 
accuracy is worst right along these common boundaries where it is ©(e^^^e'^™'/^) which in turn is largest 
just at the tips of the teeth {\pm\ = k ^og{eJf^)), where the estimate reduces to 0{e]^^). 



7. Proofs of Theorems 

7.1. Exact formulae for the potentials. Beginning with the matrix N(?i;) = N^(z(;;x,i) we may now 
express 'N{w) in terms of explicitly known quantities and the error matrix F^^{w) characterized in Proposi- 
tion 6.5 For example, if \uj\ is sufficiently small then for > we have 

Niw) = Oiw) = Piw)ta, = QMPrn'H^'^i = E± («;)Q± (u;)Pr (^^)*'Ti, (7.1) 

while for < we have the same formula just replacing P{w) and P™*(w) by — P(w) and — P$^*(w) 



respectively. But from (3.7) we see that 

sgn(5{«;})Pr (^) - (l + ^—^<J3i^wY/' + Oiw)^ (-l)"+^*a2, w ^ 0, (7.2) 
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m+l 
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- yi\ 






\y 


- yi\ 






\y 


- yi\ 





0{e]i'\y-y,\ 



-1' 



Figure 6.5. Estimates e+ (j/, s; ejv) /or expansion coefficients o/E+ (w) — I wii/i (y, s) € 
and e~^]^(?/, s; eAf) /or expansion coefficients of E^^j^(u') — I with {y,s) e ^^m+i- ^/je 
boundaries of the two interlocking regions are drawn with solid black curves. 



where 5 :— (— w*)^/^ > 0. Therefore, in terms of the expansion (6.721 of 'Ei^^{w) for small \w\, and the 
corresponding expansion of Q^(w): 

we find the following formulae for certain expansion coefficients of N(u;) (see ( 1.26| )): 

[lN%]ix,t)^[^Ei][^Qi]i~ir+Ha,, 

and 

rO-NTB j-\-irOi\.T0 u™ J-^ _ rO/S± i-l rOTn± i-l rOTpi i rOrS± i _ , _ rOrS±i-irO/S±i 



[^N%]{x,t)-'tiN%Kx,t) = a2[°Q^]-MSE±J-M?E±][°Q±]a2 + a2[°Q±]-M;Q^]a2 - ^-^a,. 
Similarly, if |u;| is sufficiently large, then we have 



and from (3.7) we have 



'(w) = I + (1 - 2m)Sa3{-w)-^'^ + 0{w'^), w ^ oo, 



so in terms of the expansion (6.721 of E^(u') for large \w\, and the corresponding expansion of Q^(w): 



(7.3) 
(7.4) 
(7.5) 

(7.6) 
(7.7) 
(7.8) 
(7.9) 



we obtain an exact formula for the coefficient [°^N^](a;,i) (see (f.25)): 

nN«^](:r,t) = nE±] + nQ±] + (1 - 2m)Sa,. 
Combining with ( |1.27p from Definition [lT] we obtain the exact formulae: 

cos(iu^(x,t)) = (-l)™[[°E±][oQ™]]i2 and sin(i^^(a;,i)) = (-l)'"[[X][oQ™]]22- (^-lO) 
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Similarly, combining (7.5) and (7.9) with (1.28) gives 



du 



N 



dt 



1 21 ■ 



(7.11) 



Now suppose t hat ( x, t) are such that (y, s) 
apply Proposition 



-2/3, 



6.5 



(x, f), s(x, t)) € with y bounded. Then we can 
and the boundedness of [qQ^] — Qm(0) ^^-^ its inverse as guaranteed by Proposi- 
tions |6]2] and [13] to obtain: 

cos(iuAr(a;,t)) = C^„^{x,t)+0{e^^{y,s;eN)) 



all holding for {y, s) € fi^j with y bounded, where 
:=(-1)"SQ^]i2 



sin(iuAr(a;,t)) = + 0{e^^{y, s;eN)) 

(x, i) = G%(x, t) + 0(e± (y, s; e^)). 



(7.12) 



Eat 



(7.13) 



1 



[[oQm] [lQm]]21 



[°lQm]l2 



[oQm]2l[?Qm]ll 



[oQm]ll[lQm]21- 



(In the last equality we used the fact that [qQ^J = Qm(0) has determinant one.) Since Q^(w) = ¥^{w) = 
G^(i(— w)^/^), the necessary expansion coefficients of Q,^(w) (see (7.3) and (7.8)) can be obtained directly 



from the partial-fractions expansions (6.37) and ( |6.58 ). In this way, we get 



implying that 







+ -J- («f^2[b+ 








[?Q™] = -^H<^2[2b+- 


5d+],2b+ -*5d+) 


(7.14) 






-i (-icr2d+,d+) , 








= (-1)" 




















= (-1)" 


-. ^m,2 *'^m,2 




(7.15) 






h 25*6^ -I- 2 *^m,2^m,l] ' 





where the components of the vectors b+ and d„j are given by (6.40) with (6.41) and we recall that S := 
(— w*)^/^ > 0. In the same way one sees that 

loQm] = I + — (a™ - «'S'c- , -icr2fc - i^c"]) 



1 



(2a„ -iSc 



-i (c,7,,icr2C,7j , 



implying that 



m+l 



HQ; 
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(7.16) 



a. 



■m,l , 'Wn,! 



(7.17) 



where a„ and c„ are given by (6.61) with (6.62) and we have S := (— w*)-^/^ > 0. 
Finally, substituting from ( 6.40 )-( 6.41) and (6.61 )-( 6.62 ) gives 



m-\-l 



2Ri 



N,mi^^t) = (-1) 



G 



N,mi^^t) = 16W* 



[R^? + 1 



(7.18) 



where 



([J±]2 + 16w;2)2 + 64?«2( j± + SKf^y 



R: 



[J±]2 + 16wl 



(7.19) 



and where J+ and K:^ are defined by (6.33) while J„ and if„ are defined by (6.54). 



7.2. Proof of Theorem |1.1| (Initial accuracy of fluxon condensates). The fluxon condensate u^ix, t) 

2 /3 

is an even function of x, so it suffices to consider x ~ Xcnt > 0. The line segment t = with Aa; = C'(£jv ) 
is mapped under {y,s) = {e~'^^^r{x,t), s{x,t)) into the union , assuming that ejsr sufficiently small. 

More specifically, it is mapped into the part of fl^ U 0,^ with y = 0(1) and po — 0{e]^^). Using the facts 
that Uo{y) — Vi{y) = 1 and Vo{y) = Wi(y) ~ — 2//6, and hence also Ho{y) = Hi{y) — ?/2/24, it is easy to 
see from the definitions (6.33) and ( |6.54 1 (using also h^{wi,) — K^{Wi,) — —'kv.jtW {w.j,)) that for t — Q and 
Aa; = ©(e^/^). 



^0+ 



while 



0{e^^) and Jf = 0(e]^/') 

Ti'- — I /n/'^i/3\ 



K+ = -2S + 0{ej^') and i^^" = -25 + 0(e^"), 
where S := (— w*)^/^ > q According to ( 7.19[ ) we then have 



R^n 



-l + 0{e]^^) and 



-l + Oie 



1/3 ^ 
N ) 



It follows from these facts and (7.18) that for such x and t ^ 0, 

C+o(x,0) = l + 0(e]v^') and C],^,{x,0) = 1 + 0{e]l^), 



S+,{x,0)^O{e]i') and S], ,ix,0) ^ 0{e'f), 



and 



G+ (x,0) = - 5 



5* 



0(e]v/=^) and G^,(x,0) = -(5 



S 



-AT 



(7.20) 
(7.21) 

(7.22) 

(7.23) 
(7.24) 

(7.25) 



Furthermore, since is an analytic function of {x,t) with (a;crit, 0) = —1, we have = — l + C'(Ax,t) = 
-1 + 0{Ax) = -1 + Oie^J^). Therefore, we may equivalently write the above formulae in the form 

G+,{x,0) = -2 + O{e]^') and Gj,,{x,0) = -2 + O{e]l'), Ax = 0{e'f). (7.26) 

Also, since the function G(-) in the initial data is differentiable at a^crit and G(a:crit) = ^2, we can also write 

G+^^{x,0) = G{x) + O{e]/^) and G]^^,{x,0) = G{x) + 0{e](^), Ax ^ 0(e^j(^). (7.27) 



Finally, recalling (7.12 ), the fact that e(j"(j/, s; eat) — 0{e]^^) for (y, s) G fig with y bounded and po — 0{e 



and the fact that (y, s; eA?) = 0{e]^^) for (y, s) e fi^ with y bounded and po — ©(e^"*), we arrive at 



l/3x 



un{x,0)^O{4'') (mod 4^) and eAr^(a;, 0) = G(a;) + ©(e^^/'), Ax = 0(e^/'), 



at 



(7.28) 



which completes the proof of the Theorem. 
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7.3. Proof of Theorem 1.2 (Main approximation theorem). Firstly, we wish to simplify the expres- 
sions for C^j^{x,t) and S^j^{x,t) given in (7.181 by localizing near the critical point Ax = t = 0, which 
allows us to write the coordinates r{x,t) and s{x,t) explicitly in terms of x and t at the cost of some small 
error terms according to (4.381. By analogy with the definition (6.11) of p,n for general to G Z, we define a 
shifted and scaled coordinate qm equivalent to t by the relation 

2m 



t = -^e^v log(e^r ) + tNqm, meZ. 



(7.29) 



Recall the spatial coordinate z exactly proportional to Ax as defined by ( 1.43 ). We will now show that there 
is essentially no additional cost (beyond that already introduced in approximating E^j(u))) in explicitly 
introducing Ax and t as coordinates. 



Proposition 7.1. Suppose that {y, s) € fl^ with y bounded and m fixed. Then 

C^N.mi^^ t) = C!N,mi^^ + C(e+ (y, s; ew)) and Sjf „^{x, t) = 5^^„,(a;, t) + 0{e+iy, s; ejv)) 
where 

Ctn(^,t) := (-1)^ 

and where 

R+ 



27?+ ■■ (7?+l2 - 1 

^m+l ^^^m S+(x,t):={-l)'^+^^^ 

_24(l-m)j^-2m/3gg„^^(-^) 



2-8m^-2(l+2m)/3g^/3g2g„(2i/„(z)ZY„(z) -Z^;,(z))2 + 16 ' 

Similarly, suppose that [y, s) G il^ with y bounded and m fixed. Then 

^N.,mi^^ t) = C'jv,m(2;, t) + C(em(2/, s; <^n)) and Sjf^^{x, t) = %,„(a;, t) + 0(e;;(?/, s; ejv)) 
w/iere 

-1 [-^m]^ - 1 



id w/ie 



2i?„ 



a'^c? Vm(x,t) := (-1)' 



[i?m]2 + 1 

_24"v2(m-l)/3g-g„-ii;^(-^) 



[R7nY + 1 



28(r„-l)^-2+4m/3g^/3g-2g„_i(27?„(z)V™(z) ^ V„,{z)Y + 16 ' 



(7.30) 
(7.31) 

(7.32) 

(7.33) 
(7.34) 

(7.35) 



Proof. First, we localize near criticality by isolating the terms in R^ depending "slowly" on x and t and 
replacing x with Xcrit and t with zero in these terms. The "slowly varying" quantities are w.^, — 'w^,{x,t) and 
the following: 



4w^,h^{w^,) 

Aw,h^{w,)W'{w,) - 4w,hf^{w,)W"{w,) - Ah^Jw*)W'{w,) 



(7.36) 



Recalling the definitions of and given by ( |6.33| ) and (6.54 1, one obtains from (7.19) the exact 
expressions 



R+ 



R^ — 



fi+eP-^Um{y) + l+^ji\^-Bm,12{y) 



\ + 1/3 



B.m,12{y)] 

[a;ne]l^e-P^-^B„,,2i{yW + ISw'*' 



where Bm,i2(2/) = 2Hm{y)Ujn{y) - l^'miv) 



/3„e-f™- V,„(y) + 7,n4^'e-P'"-B„,2i(2/) 



(7.37) 



where ^^,21(2/) = V^{y) - 2H^{y)V^{y). 
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Now, for each to G Z, a^, 7^, and are all independent of eAr, depending only on {x,t). Moreover, 
they are all analytic functions of {x,t) near criticality. Computing their limits at criticality using (4.55) and 
(|4.61|) together with (|6.18|) and (|6.42|) we obtain: 



/3± = -22±(2-4r«)^(-l±(l-2™))/3 

7± =0(Aa;,i), 



(7.38) 



and for we have the corresponding expression (4.56 1. Since by assumption to and y are bounded, it follows 

Then, since the Jacobian matrix 



immediately that s{x,t) — 0(eAr log(e^"'^)) and that r{x,t) — 0{€'Jj^) 



d(r, s) / d{x,t) is diagonal and invertible at criticality according to (4.38), it also follows that Aa; = C'(e^'^) 
and t — 0{e]\j log(e^^)), and consequently all error terms of the form 0{Ax, t) may be replaced with ©(e^"^). 
Therefore, 



Rt = 



R.„, — 



[_2-4™^-(l+2m)/3,V3gp„^^^^^(y)(^ ^ 0{e'f))]2 + 16(1 + 0{e'^')) 
[_24(™-i)z.-i+2™/34/3e-P™-ii3„,2i(2/)(l + 0(e5/'))]2 + 16(1 + 0{e'^')) 



(7.39) 



The denominators are bounded away from zero and, according to Proposition 
have simple poles at the points of 3^{Um) = >^(V„i) implying that 

ejv^'ef" 5,^,12(2/) = 0(1), for (y, s) e n+ 
e]^''e-P'-'B^My) = 0{1), for (y,s) e n^. 



5.2 



Bni,i2{y) and Bm,2iiy) 
(7.40) 



Therefore, the formulae (|7.39|) can be written in the form 



_24(l-m)j,-2m/3gp„^^(y) 



2-8m,,-2(l+2m)/3g^/3g2p„5^_^2(y)2 + 16 
_24m^2(m-l)/3g-p,„_iy^(y) 

28(m-l)j,-2+4m/3gV3g„2p„_i5^ 2l(2;)2 + 16 



(7.41) 



0(e 



2/3 
N f 



V„.(y)) + 0(e?/') 



The explicit terms we have retained in i?^ have the same form as i?^ defined by (7.32) and (7.351 but 
written in terms of the "fast" variables Pm, Pm-i, and y instead of qm, qm-i, and z. But, recalling the 
definition (7.291 we deduce from (4.38) that 



Pn 



0(e^iAa;2,e^4Ax,e^42) 



(7.42) 



and from this it follows (using ( 7.40 1 and the fact that the denominators are bounded away from zero) that 

-2^(1- 



1 — m ,,— -i 



Rl = 



R„ 



V'"/'eg'"Z^^(z/) 

2~8mj,-2(l+2m)/3g^/3g2q,„5^^^2(y)2 + 16 
_24m^2(m-l)/3g--<;™-ii;^^(j^) 



28(m-l)j,-2+4m/3g^/3g-2g„,_i5^^ 2l(2;)2 + 16 



0{e]/'eP-U„.{y)) + 0{e'^^) 
+ 0{e]i'e-P-^V^{y)) + 0{e'^'). 



(7.43) 



Also, from the definition (1.43), the definition y = r{x,t) / e^f" and (4.381, we see that 



Since the poles are simple for U„i and Vm, 

UUy)-Ura{z) = 0{eTu'^{y + 0{eT))) = 0{e%% + 0{e'^') - ^{U„,)\-^) 
V™(y) - V„(z) = 0{eTV^{y + 0{eT))) = 0{e'J,% + 0(ef ) - £?>{y^)\-^). 
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-2/3 



2/3,2n 



2/3n 



(7.44) 



(7.45) 



Now, \y - 3^(Um)\~^ = for {y, s) e n+ U fl^, so 



1 



1 



\y-^iU^)\ 



\y + Oie^) - 3^{pi^)\ \y - \y + ©(^Z^) _ ^{u,^)\ 

1 \y-^{U m)\ 

\y-,^{U^)\ \\y-^{Um)\ - 
1 1 



< 



OieT) 



(7.46) 



\y-^{llrn)\l-\Oie]/^)\ 



SO that W„(y) -U^{z) = 0{4^\y - 3^{Um)\-^) and similarly V„,(y) - V^{z) = 0(e^/'|j/ - ^(V„0|-2). For 
-^m 12 s-iid 21 the poles are double, and it follows by the same reasoning that i?m, 12(2/)^ — Bm,i2{z)'^ = 

0{e]l^\y - and B,„,2i(y)^ - B,„_2i(z)^ C'leJ/^ly - ^(V™)r^). From these considerations we 

have that 



(7.47) 



+ 0(e2,e-3p,„_i |y _ ^(v„)|-5) + o{. 



V„,(y)) + 0(e^/'), (2;,.s)ef2-, 



where i?„j is defined by (7.32) and is defined by (7.35). By comparing with the definition (6.691 of 
e+ (y, s; cn) for (y, s) e 17+ and the definition (6.70) of e~ (y, s; e^v) for (y, s) € we see that 

< = 7?± + 0(e± (y, s; ejv)), (y, s) e f7± . (7.48) 

Finally, we insert the estimate ( 7.48 1 into ^) ^^"^ "^w m(-^' ^) given by ( 7.18 ) noting that derivatives 

of the rational expressions in i?^ that appear are uniformly bounded to complete the proof of the proposition. 

□ 

Combining Proposition |7. 1| with the formulae (7.12), wc obtain the estimates 

cos{\uN{x,t)) = C^.^{x,t) + 0{e^{y,s\tN)) and sin(iwAr(a;, t)) = S"^ ^(x, t) + 0(e^(y, s; cat)) (7.49) 

holding uniformly for (y, s) £ H.^^ with y bounded. 

We are now in a position to define the multiscale asymptotic formulae for cos{^un{x, t)) and sin-{-^U]y{x, t)) 
mentioned in the statement of the Theorem: 

C{x,t;eN) := C^^mi^,^) and S{x,t;eN) := S^^mi^it), whenever {e]^^^^r{x,t), s{x,t)) e n^. 



(7.50) 



Using ( 7.49 ) we find that 



cos(iwAr(a;, t)) = C'ix, t; en) + 0(e^(e^^/^r(x, t),s{x, i); e^)) = C{x, t; en) + 0{e]^^) 
sin(i'UAr(x, t)) = S{x, t; e^) + C'(e^(e^^^^r(x, t),s{x, t); en)) = S{x, t] sn) + 0{e]^^) 



(7.51) 



hold whenever (x, t) are such that {ej^^^'^r{x, t),s{x, t)) € fi^. The exponent of 1/6 is the worst-case exponent 
as explained in the paragraph following the proof of Proposition |6.5[ On the other hand, it is clear that this 
error estimate is uniformly valid whenever m is bounded and e^y r{x, t) — y is bounded, or as is the same, 
whenever t = 0{eN log(e^^)) and Ax = 0(e^^). This completes the proof of Theorem 



1.2 



7.4. Proof of Theorem |1.3| (Superluminal kink asymptotics). Choose any integer m with |m| < B, 

2/3 

and let if > be fixed. Suppose that {x, t) lies in the rectangle TZm defined by the inequalities | Ax| < Ktjlj 

and \t — Imejv log(e^^)| < ^eat log(e^^). Because y = z + 0{e'Jj^) and Pm — Im + 0{e'Jj^), for any strictly 
positive choice of the parameter k > entering into Definition 6.1 the image I{TZrn) of TZrn in the (y, s)- 
plane under the mapping y = r(x,t) j and s — s{x, t) will be contained in the horizontal strip fJ^j U ^^+1 
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(see ( |6.14 1) with y bounded, as long as eat is sufficiently small. Also, for {x,t) mapping into I{TZ 



0+ (respectively into I(TZ„i) n fl^^^i) wc have the uniform estimate e+ (ejY^^'^r(x, i), s(x, i); ejv) = ©(e]^") 



n 



(respectively we have the uniform estimate e~_|_]^(e^^^'^r(x, t), s(a;, t); eat) — ©(e^")). 

Moreover, given an interval \ z-, z +\ with 2iy^^^\z± \ < K on which log is bounded, if the parameter 

is chosen sufficiently small, the sub-rectangle TZ^^^ C TZm defined by 



S > entering into Definition 



6.1 



\meNlog{ejy^)\ < 



the inequalities z- < z < z+ (recall the definition (1.43) of z in terms of Ax) and \t ^, 

log(e^^) will be mapped into a subregion 7(7?.^"^) of HjZm) in which both inequalities \y — 0^{Um)\ > 5 
and \y — ^^{Vm+i)] > 6 hold. It then follows that for {x,t) mapping into I{TZ'^^) n il^, the estimate 

'-1=^ n the 



0(e^ ) holds uniformly while for {x,t) mapping into I{TZ'^- 



6.5 



estimate e^_^_i{epf^^^r{x,t), s{x,t); = C(ejv'^) holds uniformly. See also Figure 

Suppose that {x, t) corresponds to a point (y, s) e £7+ . We wish to neglect the first term in the denominator 
of .R+ as defined by ( 7.32| ); we write 7?+ in the form 



R+= -^"y-'nyy- ^ ,f^O(eTe'^-B^,M)^0{eVr'e'i^\z-^^{U^)\-'), (7.52) 

where Tk is defined in terms of t, e^r, v, and |Z/^m(2;)| by (1.42). Note that since z = y + 0(e5/^) and 
9m = Pm + C'(ejv'^)' ^^'^ since \y — ^{Um)\ > ^^jv^ for {y, s) G 17+ , the upper bound for / m ay b e replaced 
with / 



gFm|y _ ^(U^)\-^). Therefore, as pointed out in the proof of Proposition 7.1 / = 0(1) 
uniformly for {y,s) E $1+, although it is typically smaller pointwise. By using this result in (7.31) and 
recalling (7.49), it follows that whenever {x,t) E TZ,n maps into 11+, 

cos(iujv(^,<)) = (-l)'"sgn(Z^,„(z))sech(TK)+0(e-l^-l/2) + 0(e+(e-'/'r(x,t),s(x,i);eA,)) 



sin(iujv(x, t)) = (-1)™+! tanh(rK) + Oie-^^^-^f) + 0(e+ (e^'/3r(x, t),s{x, t); e^)) 



(7.53) 



As shown at the beginning of the proof, the error terms involving e+ are either of mag nitude Oie]^^) 
and in particular o(l), or, if {x,t) E 7?.™^, of smaller magnitude 0{e]^^). The remaining error terms 
are bounded, and they are negligible unless both Tk = and / fai ls to be small. But since / 

"^(Z^m)|~^) and {y,s) E f2+, we see from Definition 



He 



1/3, 
N 



0{e]^^), or equivalently, e 



-2m/3 t/e 



0{e 



1/3n 
N ) 



6.1 



that / fails to be small only if 
Now, under this condition we have from 



l/3^ 



N 



( [02] ) that Tk > 0(1) + i log(e^^) -flog |W™(z)|, so to have Tk = 0(1) requires that \z-^{Urr, 
Moreover, if {x,t) E and its image (y, s) = {e~^'^r{x,t),s{x,t)) lies in 17+, e"l^'<l = 0(e-l9'"l) and 

/ = 0(e]^^e'''"). The inequality \q^\ < | log(e^"^) then implies that if {x,t) E Tl'^^ with image in 17+, then 
g-|Tk|j2 g-2|TK|y2 \jQfj^ dominated by the 0(e]^^) estimate of the error terms involving e+. 

Now suppose instead that {x,t) corresponds to a point {y,s) E 17^+1. As before, we wish to neglect the 
first term in the denominator of R^-^^i as defined by ( 7.35[ ), so we write ^,^+1 in the form (using the identity 



R 



m+l 



sgn(Z^m(2:))e '^^ 1/3 _„ , , 1/3 , s.-ix ,„ 
f^2 ' " = C'(ej^ e ''"B™+i^2i(2:)) = 0(e^ e '^'"\z - £!f(Urn)\ )■ (7. 



54) 



Again, the error bound for h may be replaced with 0{e](^e ^"ly — 3f{Um)\ ^) and so h — 0(1) holds 
uniformly for (y, s) E ^l^^^; it then follows by using this result in (7.34) and recalling (7.49), that whenever 
{x,t) E TZm maps into 17^^^^, 

cos{^UNix,t)) = (-l)™sgn(Z^„(z))sech(rK) + 0(e-l^^l/i2) + 0(e,-^^(e^2/3^(x,i),s(x,i);ejv)) 

sin(i7.^(a;, t)) - (-1)'"+^ tanh(TK) -t- 0{e'^^^>'^h^) + 0(e„+i(e^'/V(x, t), s{x, t); ejv)). 

In particular, we note that the leading terms are identical with those arising from the analysis in the case 
that the image of {x^t) lies in 17+. By completely analogous reasoning as in that case, one learns that the 
sum of error terms is o(l) unless both Tk = and also h fails to be small, which occurs in this case only 



(7.55) 



if eP" 



1 /3 

0(ej^ ), or equivalently, if e'' 



2rn/3 t/eM 



1 /3 ^ 

0(e^ ). Under this condition, the definition (1.42) 
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implies that Tk < 0(1) - 5 log(e^^) +log \U„i{z)\, so to have Tk = 0(1) requires that \z-^{Um)\ = C(e]v^^)- 
Moreover, if {x,t) e Tl^^ with image in then both error terms are ©(e]^^). This completes the proof 

of Theorem 11.31 



7.5. Proof of Theorem 1.4 (Grazing kink collisional asymptotics). Recall that zq £ ^{tlm-i) is a 
simple zero oilAm-i, and hence also a simple zero of Vm+i- Note that from the recurrence relation (5.62) it 
foUows that V^_|_i(2;o) ~ —l/^m-iC^o)- Since z = y + 0{e'^^'^) and t ^ s + 0{e^^) for the range of x and t 
under consideration, the region in the (a;, i)-plane defined by the inequalities 

|t-(|m-|)eArlog(e^^)| < ieArlog(e^^) and |z-zo| < //e]^/^ exp(|t- (|m- |)eAr log(e^^)/(2eAr)) (7.56) 

is, for both /i > and e^r sufficiently small, mapped into the union ri^j_]^ U ^^+1 ^^e (y, s)-plane, where 
{y, s) = {ej^^'^r{x, t), s(x, t)). In particular, each point of the image lies either in the upward-pointing "tooth" 
of i^m-i centered at y — zq or in the downward-pointing "tooth" of centered a.t y — zq. It follows 

that for the points whose images lie in the error term e^_i{ej^j^^'^r{x,t), s{x,t); e) is 0{e]{^) while 

for the points whose images lie in the error term e~^]^(ejy^^^r(x, t), s{x, t); e^r) is also 0(e]v/^) so that 

in particular these error terms arising in (7.491 will always tend to zero with ejv- Moreover, the inequality 



1 /3 

|z — zqI = 0{^N ) implies that \y 



zq\ 



1/3 

0{ej!j ), from which it follows that under this condition {x,t) will 



be mapped into the parts of and ^^^j^i (the centers of two matching teeth) in which the better error 

estimate of 0{e'J,^) holds for e^_i and e~_|_;^ respectively. 

Suppose that the image of (a;, t) in the (j/, s)-plane lies in the upward-pointing tooth of centered at 

y = Zq. In this part of the image, the uniform estimate e^*^ = ^(1) holds. By Taylor expansion. 



U„.,,{z) - U'^-i{zo){z zo){l + 0{z zo)) = %^X»-i(^o) ^^^^ ^ _ ^^^^ 
B,n-iMz) = 2H„,-i{z)U^_i{z) - U'^_^{z) = -U'^_^{zq) + 0{z - zq), 
where the spatial coordinate Xc, is defined by ( |1.49 ). It follows from (7.32) using e^"^ = £>(!) that 

= -sgn(W;,_i(zo))^G sech(rG)(l + 0(z-zo)), 



and then that 



_,{x^t) = (-l)'"-^sgn(^;,„,(zo)) , +0{z-Zq) 

1 + sech (Tg) 



(7.57) 



(7.58) 



(7.59) 



+ O{z-zo). 



H-X2sech2(rG) 

Combining this formula with the above estimates of e^_i{x, t; ejv) and the formulae ( 7.49 ) proves the Theo- 
rem for those {x,t) whose image in the (y, s)-plane lies within the upward-pointing tooth of ^^^^ centered 
at y = Zq. 

Now suppose instead that the image of (x, t) in the (y, s)-plane lies in the downward-pointing tooth of 
^m+i centered at y = zq. In this part of the image, the uniform estimate e""^*^ — 0(1) holds. By Taylor 
expansion and the fact that V4i-|_i(zo) = ^l/^m-iC-^o), 



Vm+l{z) 



1 



Ka-liz^) 



(Z-Zo)(l+O(z-Zo)) 



^1/3 



2vum'„,_,{zQ) 



Xg(1 + 0(z-zo)) 



(7.60) 



Bm+iai{z) = - 2iJ„i+l(^)Vm+i(z) 



^m-l(^o) 



+ O{z-zo). 



It then follows from (7.35) using e ^'^ = 0(1) that 



R;^^^ = sgn(Z^;„_i(zo))XG sech(rG)(l + 0{z - zo)), 
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(7.61) 



and then that 



CM,m+ii^:t) - (-l)"'-^sgn(^:„_i(zo)) ^ l^^r^Sfl ^ +g(^-^o) 

1 + scch (Tg) 



(7.62) 



Combining this formula with the above estimates of e^J^-^^{x, t; e^) and the formulae (7.49) proves the Theo- 
rem for those (a;, t) whose image in the (y, s)-plane lies within the downward-pointing tooth of ^^^i centered 
at y = zq. 
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